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GILilLi^ PiJPI RSDDY 
Ph.D. 

department of Mechanical Enginfeering 
Indian Institute of leohnology Kanpur 
December 1975 

AU!K)MaIjin> OPTIMUM DESIGN iiHD EEUiiBILIII AUiJDYSIS OF 
MACHINE TOOL STHICTDEES 

^ eomputational capability for the automated optimum design 
of complex machine tool structures to satisiy static jigidity» 
natural frequency and regenerative chatter stability requirements 
is developed in the present work, Mare specifically, the minimum 
wei^t design of Warren type lathe bed and horizontal knee type 
milling machine structures using finite element idealization is 
considerdd in this exploratory study. The Warren type lathe bed 
is cptimized to satisfy torsional rigidity and natural frequency 
requirements, whereas, the milling machine structure is optiiaized 
with static rigidity of the cutter centre, natural frequency and 
regenerative chatter stability •onstraints. 

A survey of the available literature indicates that the 
use of computers in the machine tool manufacturing industry has 
up to the present time been confined to the 11153 lementation of 
finite element technique for the purpose of static and dynamic 
analysis only and the potentialities of optimization techniquss . 
have not yet been exploited by the machine tool structure! designers. 
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In designing these structures, . the most unfavourable cutting 
conditions and sizes of work piece are taken into account. The 
torsional deflections and the first two natural frequencies of 
vibration are the behavior constraints in the design of lathe beds. 

In the design of horizontal milling machine structures, the maximum 
deflection of cutter centre in any direction, the natural frequencies 
of vibration and the minimum negative inphase cross-receptances of 
cutter centre relative to the table behavior 

constraints. The constraint on is oregenerative chatter 
stability of the milling- machine. 

In the design of lathe bed, the thickness, of the main 
members, thickness of the lacing diagonals, the width and the 
thickness of the flanges (stiffeners) on the main members and 
lacing diagonals and the width of the bed are considered as, design 
variables. In the case of milling machine, the thickness of over— arm, 
the thickness of column and table, tie depth of the column at the 
base, the depth of the column at the top end, the width of the 
machine, and the cross sectional area of ribs on the over-aim and 
its joint with the column are taken as design variables. 

The idealisation using triangular plate elements with a 
3 -term in-plane and a 9-term transverse displacement model and 
frame elements has been found to be effici.ent. - The triangular 
plate elements are used to idealise the main members and lacing 
diagonals of lathe feed, and ■:he colton, oirer-aim and table of 
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aiilliAg laaciiine structures. Ihe frame elements are used to 
idealise the carriage guides (guides on the lathe bed), stiffeners 
on the main members and the lacing diagonals of lathe bed. In the 
case of milling machine, the ftame elements are used to model the 
ribs on the over-arm, the oveor-arm joint with the column, the 
arbor and the arbor support. Ihe grouping of identical and 
identically disposed finite elements (having the same transformation 
matrix from elemental to global coordinate system) has resulted 
in savings of about 20fc to 80fi> of the computer time in generating 
the global stiffness and mass matrices. 

Ihe cholesky decomposition of symmetric banded matrices has 
efficiently been applied in solving the equilibrium equations and 
in obtaining a partial solution of the eigenvalue problem by 
Eayleigh-Ritz sub-space iteration algorithm. This algorithm solves 
the eigenvalue problem directly without a transformation to the 
standard form. This method is the most efficient one for solving 
a partial eigen solution of large structural systems. The 
reoeptances of the cutter centre relative to the table of horizontal 
millixig machines have been obtained by using modal coordinates 
taking the damping matrix as a linear combination of the stiffness 
and mass matrices. 

Hie constrained optimum design problem is cast as a non- 
linear mathematical programming problem. The interior penalty 
function method, with a variable metric unconstrained minimization 
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technique, is used to solve the minimum weight design problems. 

The ^somputation of the gradient and the slope of the <)-function 
has been carried out by a finite difference method which incorporates 
the rates of change of the response quantities with respect to the 
design vaxiables. The computational experience shows that the 
approximate methods used in the reanalysis cycle of optimisation 
have baen quite efficient and reliable. 


In the design of horizontal milling machine structures^ 
tho value of the minimum negative inphase cross receptance of 
the cutter centre relative to the table bas been computed 

for the most common range of cutting conditions. 


(The coupling constants for plain milling operation (for the 
most common cutting range) have been computed both by using tho 
basic theory of metal cutting and by using an empirical fomula, sug- 
gested by Vulf , which computes the average tangential force in 
plain milling. A close agreement is found between the values 
obtained from the two methods. The critical value of the coupling 
constant (r^ ^-p... ^) has been selected from this analysis to obtain a 


bound on &, 


MUT* 


The design examples of lathe bed indicate that the thicknesses 
of both tbe main members and lacing diagonals assume their minimum 
permissihle values at tne optimum point. The widths and thicknesses 
of the flaxiges on the main members and lacing diagonals also 



decreased ftoia the starting design, Ihis establishes that in this 
type of structures > thin walled structures are preferable. Ihe width 
of the lathe bed increased such that the correspondir^ lacing angle 
becane 53»5® at the optimum. This value of lacing angle at optimum 
agrees very well with the results obtained by Badauri, Moshin and 
Ihornly sdio analysed and tested a number of 7/arren type beams under 
static torsional loads. 

The optimization results of horizontal milling machine show that 
the thickness of the over-arm, the column and the table decreased ag 
the optimization progressed, Iherefore the results indicate that 
even in milling machines, thin wailed structures are preferable, Ihe 
cross sectional dimensions of ribs on the over-arm and its ;joint with 
the column also decreased. At the optimum, the over-all dimensions 
of the milling machine are: depth at the base of the column is 50.12 cm. 
depth at the top of the column is 41.5 cm, and the width of the machine 
is 28.8 cm. In the design of horizontal milling machines these over- 
all dimensions are very important. 

At initial and optimum designs, the reliability study of 

milling machine structures is made against the three response quantities 

i.e., deflection of cutter centre (d ), first natural frequency of 

c 

vibration (w^), and average height of the table (h), ijhe 

position of the cutxer centre (c), modal damping factors )»Yomg*s 
modulus and the load on the machine (p) are taken as probabiliS'fci<? 
quantities. A comparison of reliability results at initial and optimum 
designs is made. A sensitivity of the reliabilities with respect to 
the coefficient of variation of the various random design parameters 
is also conducted. A method of including reliabilily analysis in the 
optimum design procedure is also indicated. 
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It is s ^iLSi? omsxy IrO I^slss ibs issx^n of mso^ioo 

tool primarily upon tke requir«nents of static rigidity and minimum 
natural frequency of vibration. The effects of different machining 
parameters like cutting speed, feed and depth of cut as well as the 
size of the work piece also have to be kept in mind by a machine tool- 
structural designer. Eor a tentative design-, the machine tool is 
analysed for various natural frequencies, dynamic rigidily and chatter 
stability. Based on the results of this analysis, suitable modi- 
fications are made, by a process of trial and error, to satisfy the 
design requirements. Ihis procedure of trial and error is adopted 
mainly because of the complex nature of machine tool structures and 
also because of the lack of a suitable design procedure ’ that can 
handle all the requirements simultaneously. Moreover, the modifications 
to satisfy the design requirements are often not veiy extensive, Ihis 
is mainly because of the fact that the designers seek minor improvements 
over the existing designs. On the other hand, if alrogether new 
design specifications have to be satisfied, these trial and error 
methods are not useful and only unified procedures like the autoioated 
optimum design methods are a^iplicable in such situations. Hence it 
becomes necessaiy to develop design procedures that can consider all 
the design requirements simultaneously; 
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!Ehe dynamic stiffness of any machine tool structure has an 
important influence upon both surface finish and metal removal rate# 
Ihus to achieve a high quality surface finish, forced, vibration 
amplitudes must be controlled whilst in the iiiterests of high metal 
removal rate the machine should have high dynamic stiffness to 
resist the onset of self excited vibrations, commonly known as 
chatter. 


The aim of the present work is to develop a capability for 
the automated optimum design of machine tool structures to satisfy 
static rigidity, natural frequencies and chatter stability requirements 
Since the cutting conditions vaiy substantially in a machine tool, 
the necessity of analysing its structural reliability becomes important 
Hence the reliability analysis of machine tool structures is also 
considered in the present work. 

1.1 REVIEW OE HTERAKJEE 

Since the present work deals with the machine tool structural 
analysis, structural optimization and reliability analysis, a brief 
review of the available literavare in each of these fields is given 
in the following sections. 

1.1.1 Machine Tool Structural Analysis 

The use of computers in the machine tool manufacturing 
industry has up to the present time been confined to the imple- 
mentation of finite element technique for the purpose of static 
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and dynamic anadysis only. Hie first publication describing a 

foimalised digital computer method for calculating the deformation 

of machine tool structures appeared in 1964 \ In this paper, 

Taylor and Tobias described the application of a finite element 

program involving the use of slender beams to represent the 

structural parts of a radial aim drill and a lathe. Since then a 

number of such applications have appeared in the literature. ^ 

2 

idealizing the structure by frame elements, Cowley and Fawcett 
analysed a plane milling machine for static deflections, natural 
frequencies and mode shapes. The authors have studied the effect 
of flexibilities between joints on natural frequencies and mode 
shapes. They found that for those mode shapes for which substantial 
deformation between joints is involved, the effect of joint flexibility 

is maximum. The static analysis of Warren type lathe beds was made 

3 

by Badauri, Moshin and Thornley . The authors have studied, both 
analytically and experimentally, the effects of breadth-to-depth 
ratios and lacing angles. They investigated the possibility of 
obtaining optimum stiffness-to-wei^t ratios of the Warren beams 
from their results. In reference 4, Thornley and Howes studied the 
static and dynamic behavior of Warren type beams by experimental- 
methods. 

As a co-operative programme during the period January 1970 to 
March 1971, a model milling machine was analysed for deflections, 
natural frequencies and modal shapes by a G.I.E.P. group (international 
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Institution for Production Engineering Research) in which professors 

from several European universities and one Japanese university 
5 6 

participated ’ , They have used various types of finite elements 
like frame, triangular, rectangular and prismatic elements in 
analysing a milling machine structure. Ihe primary objective of 
the collaboration was to compare the results obtained by using 
different finite element idealizations. 

7 

In the year i960, Koenigsberger and Said concluded that 

the dfinger of self excited vibrations would appear to be slight 

in the milling machine, but resonance vibrations due to fundamental 

frequency of cutting or any of its harmonics (forced vibration) 

would have to be carefully watched. Two years later, in 1962, 

8 

working on a very similar type of machine, Andrews and Tobias 
have come to precisely the opposite conclusion, namely, that self 
excited vibration was a significant limitation in horizontal milling, 
but forced vibration was relatively unimportant. 'The authors 
distinguished between forced vibration and chatter during horizontal 
milling. !Ihey found that forced vibration amplitude is proportional 
to tht. force fluctuation amplitude, while the onset of chatter 
depends on the force variation for a given amplitude of surface 
wave . 


In reference 9, Taylor described a technique for predicting 
from design drawings the chatter stability of machine tool structures 
by computing responses to excitation from computed modes of 
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vibration and assumed damping constants. Out of several theories 

10 . 11,12 

available for relating chatter stability to response loci , 

the author used regenerative chatter theory with penetration 'rate 
effects neglected”'^. This process is illustrated with reference 
to a lathe model and three versions of a milling machine. The 
authors established the superiority of machines with box type over- 
aim over those with bar type over-aim, as characterized by hi^er 
values of minimum stability under all cutting conditions. 

Tlusty and Polacek, in reference 13, analysed the role of 
vibration of machine tool strictures on the process of chatter with 
a scope to increase their stabilities for all possible cutting 
conditions. With reference j horizontal milling machines, these 
authors recognized the "weak links" of the structure and recommended 
maximum rigidity for the over— arm and its j oint with the coloumn* 
They also recommended that the natural frequency of vibration 
corresponding to the vertical mode should be slightly hi^er than 
the one coixesponding to the horizontal mode. The satisfaction 
of this requirement ensures a good stability in down-milling because 
of the advantageous interaction of both the modes. The effect of 

certain design features of horizontal milling machines was also 

14 

investigated by Said . Koenigsberger and Tlusty gave a state of 

... • 

art discussion of tho machixie tool structural design in \3(i • 

16 

Yoshimura and Hoshi deyeloped ooaputer program for the 
computation of stiuctural dynamics -and dynamic optimization of 
machine tool structuraesf Based on the .energy equations of a 
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mechanic.'.il system at resonance, the authors developed an approach 
to the optimum design of a planar-type milling machine. 

1.1.2 Structural Optimization 

During the last few years, several optimization investigations 

have been reported in the fields of civil engineering and aircraft 

structural design. Using the method of feasible directions, Pox 
17 

cind Kapoor reported a capability for the minimum weignt optimum 

design of planar truss-frame structures with inequality constraints 

on the maximum dynamic displacement, stress and natural frequencies. 

18 

iicOart, Hang and Streeter developed a steepest-descent boundary 

value method for the design of structures with constraints on 

strength and natural frequency. In reference 19, Turner attempted 

to minimize the mass of a structure with a specific natural frequency. 
20 

Zarghamee presented a method for the minimization of the lowest 

natural frequency of a structural system of given wei^t by varying 

the stiffness of its component elements. The author used the 

21 

equations developed by Pox and Kapoor , to express the rates of 
changes of frequency with respect to the design parameters. In 
reference 22, Zar^amee considered the problem of the minimum weight 
design of a planar-truss structure with a general stability constraint. 
In this work, the gradient projection method of Hosen, in conjunction 
with the steepest-descent and alternate steps method of Sebmit, was 
used for the optimum design of rhe descretized structure. 
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In the field of aircraft structural design, Schmit and 

23 24 25 

Thornton , Giles and Stroud, Dexter and Stein attempted to 

solve various optimization problems. The automated minimum 

26 

weight design of supersonic aiioraft wings was attempted by Rao 

by placing limitations on the natural frequencies, flutter mach 

number, the elastic deflection at the tip of the wing and the 

stresses induced in the skin at the root of the wing. In reference 

27, Parshi and Schmit presented a method for liie minirntm weight 

design of stress and size limited trusses under multiple load 

conditions. When applied to a limited class of structures and a 

limited number of failure modes, the authors have achieved an 

acceptaole efficiency with an opportuirLty to find a global minimum. 

28 

darasingam and Sridhar Rao obtained a minimum weigh't design of 
linear elastic cantilever and simply supported trusses by solving a 
sequence of convex prograiaming problems. The authors used linear 
constraints involving a smaller number of design variables instead 
of solving one math ana tical programming problem with non linear 
constraints in a larger number of design variables. In reference 
29, Twisdale and Kachaturian presented a generalized dynamic 
programming approach for the optimization of planar structural 
systems. 

1 . 1.3 Reliability iinalysis of Structures 

When the parameters affecting the, strength of a structure 
and the loads acting on it are statistical in nature, the conventional 
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design approaches based on the concept of "factor of safety" cannot 
be used to maintain a proper degree of safety. A more rational 
criterion, in the presence of random design parameters, will be to 
base the structural design on the concept of reliability or 
probability of failure. Reliability analyses in structural 
engineering recognize that both loads and strengths have statistical 
frequency distributions that must be considered in evaluating safety. 
Since the design parameters like cutting conditions, dimensions of 
the work piece and the location of the tool are random in nature in 
machine tools, an analysis, based on the principles of reliability, 
becomes important. 

30 

Rreudenthal ,.in the year 1956, explained that the most 

rational way of describing the overall safety of structures is in 

terms of reliability or probability of failure. In reference 31 y 

Moses and Rinser have demonstrated that an overall level of structural 

safety can be prescribed in teims of a rational criterion like 

probability of failure, and minimum weight structures can be designed 

to meet the prescribed safety level. In the year 197^, Moses and 
32 

Stevenson considered the subject of sensitivity of statistical 

parameters and presented methods of incorporating reliability 

analysis into optimum design of trusses and frames. Since then a 

33 34 35 

nuimber of such applications have appeared in literature ^ ^ 

which show that a probabilistic design is a practical possibility. 
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Reoeiatly, the probabilistic design conoeptb have also been 

36 

applied in the design of mechanical systems. Mischke presented 
a formal relationship between the reliability and the factor of 
safety of a mechanical element. Ry treating the factor of safety 
as a random variable, he used Bienayme-Ohabyshev and Gamp-Meidell 
theorems to derive expressions for the mean and the variance of the 
factor of safety for any specified value of reliability. In 
reference 37, Rao has developed a probability based design method 
for the desigi of mechanical power transmission systems like gear 
Trains. By idealizing the transmission system as a weakest-link 
kinematic chain (similar to a weakest-link structure), the design 
has been made to achieve a specified reliability with respect to 
bending and surface wear modes of failure. 

1.2 OBJECTIVE AND SCOPE OE THE PEESEHT WORK 

It can be seen from the available literature that the poten- 
tialities of optimization techniques have not yet been exploited 
by the machine tool structural designers. This mi^t be because of 
the complex nature of machine tool structures and the troublesome 
behavior constraints that arise due to the wide range of operating 
conditions in machine tools. In the present woik, the automated 
minimum weight design of Warren type lathe bed and horizontal knee- 
type milling machine structures is considered. In designing these 
structures, the most unfavourable cutting conditions and sizes of the 


10 


work piface arts takea into account. Dae torsional defl<^otions and 
the natural frequencies are the behavior constraints in the design 
of lathe beds. In the case of design of lailling machine stiuctures, 
limitations are placed on the maximum cutter centre deflection, 
the natural frequencies of vibration and the minimum negative inphase 
cross receptance of the cutter centre relative to the table. 

In order to predict the static and dynamic characteristics 
of a machine tool structure accurately and consistently, the 
structure has to be represented by a suitable model. Ihe model 
raust be simple enough from computational point of view for embedment 
into the iterative design-analysis loop of the optimization procedure. 

In the present work, the finite element displacement method 
is used to model the machine tool structures. There are several 
ways of modelling the structure by the finite element method. There- 
fore it becomes necessary to make a comparative study of the 
various finite element idealizations to determine a suitable model 
for use iti the optimization work. In this work, a study is under- 
taken to find the relative efficiencies of two different fiiiite 
clement idealizations for modelling machine tool stiuctures. The 

idealization using triangular pj.ate elements with a 5 -term in-plane 

38 39 40 

and a 9 -term transverse displacement model ’ ^ and frame elements 

has been found to be efficient. The triangular plate elements are 
used to idealize the main members and lacing diagonals of lathe beds, 
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and the colouiaa, over-arm and table of milling machine structures. 

The frame eloraents are used to idealize the . carriage guides 
stiffeners on the main members and on the lacing diagonals of lathe 
beds. In the case of milling machines, the frame elements are 
used to model the ribs on the over-arm, the over-arm joint with 
the coloumn, the arbor and the arbor support. In the present work, 
identically oriented finite elements of the same size and shape 
(having the same transformation matrix between local and global 
coordinate systems) are grouped together in generating element 
stiffness and mass matrices. This resulted in savings of about 
20 to 80 percent of the computer tine in generating the global 
stiffness and mass matrices. 

Ibr the dynamic optimization of large structures using finite 
element method, a designer is generally confronted with two probltnis, 
namely, the computer storage and the computer time. The eigen value 
problem h:;is been solved by using one of the most efficient solution 

techniques developed by Bathe and Wilson for large structural 

41 

systems . In this technique the Hayleigh-Ritz sub-space iteration 
algorithm, which solve;^ the eigen value problem directly without a 
transformation to the standard form, has been used. In this work, 
the Gholesky decomposition of symmetric banded matrices, storing 
only the upper triangular matrix, is used for solving the equilibrium 
equations. By using a judicious descretization and node numbering 
scheme, it has been possible to reduce the band width of the stiffness 
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matrix of the structure. A snaller baud width, apart from reducing 
the computer storage, considerably reduces the computer time for 
static and eigen solutions. 

Hie frequency response of the structure has been obtained 
by using nodal coordinates, and by taking the damping matrix 
proportional to a linear combination of stiffness and mass matrices. 
Since the present day knowledge is not sufficient to estimate the 
modal damping factors from the blue prints of a given structure, 
the values of damping factors have to be obtained from experinental 
results on similar structures. In the present work, equal modal 
duraping factors of value 0.06 have been used for the first few modes. 

In machine tools, chatter occurs due to the interation of the 
cutting forces and the machine tool structural dynamics. To derive 
the limits of chatter stability for the machine tool and the cutting 
system, it is necessaiy to assume a force relationship for the 
cutting process and relate this to the machine dynaiiics. In this 
woik, a simple dynamic cutting force relationship, assuming a direct 
proportionality between the force and the undefoitaed chip thickness, 
is taken. Fo definite criterion has been established so far 
for taking the critical proportionality constants 'between the force 
emd the undeformed chip thickness in plain milling operation. In 
the present work, a study of proportionality constants in plain 
milling has been .lade to choose a critical value for incorporating 
in the chatter stability constraint. 1 metnod of including the 
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rainiwum Uiigattv© io^hase qross-receptance of cutter centre relative 
to table (in tlae design) is also developed in this work. 

Ihe constrained optimua design problem is cast as a nonlinear 
mathematical programming problem. The interior penalty function 
method, with a variable metric uncons t3rained minimiaiation technique, 
is used to solve the minimum weight design problem. The ijdnimizing 
step lengths in the unconstrained minimization are detenained by- the 
cubic interpolation method. The computation of the gradient and 
the slope of ^-function has been carried out by a finite-difference 
method which incorporates the rates of changes of the response 
quantities with respect to the design variables. In the design of 
the lathe bed, the main dimensions .of the lathe bed, the thickness 
of the main members and lacing diagonals, and the cross sectional 
areas of thoir stiffeners are taken as design variables. In the 
case of milling machine, the ot, orall dimensions and thickness of the 
coloumn and the over-ara, and the cross sectional dimensions of 
ribs on the over-arm and on its joint with the oolouim are assumed 
to be the design variables. The optimum design of a lathe bed 
and a milling machine structure are considered to' illustrate the 
effectiveness of the method proposed. 

In the optimum design of machine tool structures, all the 
design parameters have been assumed to be deterministic* However., 
in actual practice, most of the design parameter^ lice the cutting 
conditions, the material and the structural properties vary considerably' 
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Hence the deyelop/aeat of a more relastic analysis procedure, which 
t'Oices the variability of the design parameters into account, becomes 
important, Ihe present work aims at developing a method for the 
reliability analysis of machine tool structures by considering the 
various design paraaueters as random variables with known probability 
distributions. The reliability against any specified criterion 
like static rigidity, frequencies and dynamic' response of the machine 
tool structures has been analysed by treating the machine tool 
structure as a single member acted on by a single load. All the 
random variables are assumed to be normally distributed for simplicity. 
The linearization using Taylor’s series expansion (about the mean 
values of the random variables) has been used in deriving the 
expected value and the variance of a function of several random 
variables. The reliability analysis of the milling machine 
structure, considered in the optimization work, has been presented 
to illustrate the method developed. 

1.3 ORGAI'TIZATIOH OP THE THESIS 

The thesis is organized into 9 chapters and 2 appendices^ 

The formulation of the optimization problem is given in Chapter 2, 
Chapter 3 deals with a comparative study of two different finite 
element idealizations for modelling differ^ent machine tool 
structures. The static and eigen value results obtained for 
different' structures have been compared with those available in 
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the literature. OIhe idealization, for which the results conpared 
closely with the known results, has been finally selected for use 
in the optimization program , 

Ihe expressions for the cutting forces developed in lathe 
;,ind horizontal milling machines are derived in Chapter 4, This 
chapter gives a discussion of the regenerative chatter stability 
in horizontal milling machines* A method of including the 
regenerative chatter stability constraint in the optimization problem 
has also been developed in this chapter. Chapter 5 deals with the 
static and the dynamic analysis of machine tool structures. Ihe 
algorithm used for the eigen value analysis is discussed in some 
detail. 


Chapter 6 gives a discussion of the optimization algorithm 
used in the present work. The computation of the gradients of 
the response quantities along vdth the results of the approximate 
analyses are presented in this chapter. Some illustrative examples 
are given in Chapter 7 to demonstrate the effectiveness of the 
present approach for the design of machine tool structures under 
static and dynamic constraints. 

The reliability analysis of machine tool structures is the 
subject matter of Chapter 8, A method of incorporating the 
reliability constraint in the optimum design of machine tool 
structuies is indicated in this chapter. Tinally, the conclusions 
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drawn from the present study and the recommendations for future 
research are givwu in Chapter 9. 

Hhe generation of the stiffness and mass matrices of a 
triangular plate bending element in the global coordinate system, 
by considering both in plane and bending actions, is described 
in appendix A. Appendix B contains the description and other 
details of computer program. 



CHiiPTER 2 


POMJLATION OF lES DESIGN PROBLEM 

yy'en a means for predicting the oehavior cf any design is 
ayailable, when limitations on the performance and other external 
constraints on the design can he stated and when some acceptance 
critaria can he established, a design modific '-tion prohlem co-n he 
cast as a mathematical programming prohlem. 

2,1 DESIGN PHILOSOH-IY 

Any structural design prohlem can he formulated and solved 
according to either deterministic or prohahilistic design philosophy 
If all the quantities affecting the design problem arc deterministic 
tho design prohlem can be solved according to the deterministic 
design philosophy. On the other hand, if some of tho design para- 
meters are random in nature, the design problem has to he solved 
according to the probabilistic design philosophy. 

In the deterministic design philosophy, a general mathe- 
matical programming problem can he stated as follows s 

Minimize a multivariable function f(x) subject to the constraints 

gj(X) < 0, j = 1,in 

where X is a n-dimensional x’-ector consisting 
of the variables X^, X 2 ,...,X^. 

In the probabilistic design philosophy, a general 
mathematical programming problem can be stated as follov;-s ; 
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Minimize a raultivaria'ble function f(X,Y) 

Su’bject to the prohabilistic constraints 
P ^ 1 i P j » ^ ~ 1,2, ...jm. 

•where X is the vector of design variables, and T is the vector 

of other parameters affecting the design problem. Here the 

components of X and Y are assumed to be random variables. In 

equation 2.2, f represents the mean value of the objective function, 

and P [g-(X.,Y) ^0] ^ p- denotes that the j constraint has to 
3 3 

be satisfied with a probability of greater than or equal to some 
specified quantity, p^ (j = 1,2,...,m), -where 0 1. P j * 

In the present work, only the deterministic design 
philosophy is used for the optimization of machine tool structures. 

2.2 OBJECTIVE ITJECTION 

The function f(X) in equation 2.1 is called the objective 
or criterion function and its choice is governed by the nature of 
the problem. The objective of any structural optimization problem 
is generally taken as minimization of weight or cost. . The 
minimization of weight is more frequently taken as the objective 
since weight can be easily quantified and it reflects the material 
costs and, in some cases, the manufacturing and the transportation 
costs. The minimization of the weight of the struct'ure is taken 
as the objective in the design of Warren type lathe bed and 
horizontal knee-type milling machine structures. The maximization 
of static rigidity or natural frequency of vibration or dynamic 
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rigidity of the rriachine tool structure can also he taken as 
ohjective function in some cases. 

Although, only the deterministic design philosophy is 
used for the optimization of machine tool structures, an attempt 
is made to analyse the behavior of milling machine structures hy 
taking the following quantities as probabilistic ; 

i) the cutting conditions, ii) the material properties of the 
structure iii) the modal damping factors, iv) the position of 
the cutter and v) the position of the table. Suggestions are 
given for incorporating the probabilistic analysis for designing 
machine tool stiructures according to the probabilistic design 
philosophy, 

2.5 DESIGN COlISTRAiroS 

The design requirements to be satisfied in the case of 
Warren— type lathe bed are i) the torsional rigidity must be 
greater tlian a specified quantity and ii) the natural frequenciss 
of the structure are to be excluded from certain bands. In the 
case of milling machine, the design constraints are i) the maximum 
deflection of the centre of the cutter in any direction should 
not exceed a certain prescribed value ii) The natural frequencies 
of the structure are to be excluded from certain bands, and 
iii) The machine should not chatter under the stated cutting 

These constraints are stipulated so as to achieve a 


conditions. 
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high quality surface finish, to avoid mild harmonic forcing that 
might ca.usG resonance and to increase the metal removal rate vfhioh 
may he affected due to the onset of chatter. In the design of 
metal cutting machine tools, the static and dynamic stiffness 
requirements are often more important than the load carrying 
capacity requirements since the induced stresses corresponding 
to the perinissible deformations are generally far less in value 
than those permissible for the various materials. This charac- 
teristic has been observed in the present work also. Hence the 
strength was not considered as a design constraint in the present 
work. 

2 .4 lESIGF V1EI131ES 

Once the objective function and the design constraints 
are specified, the statement of the optimization problem will be 
complete as soon as the design variables are identified . 

2.4>t bathe Bed 

The details of a Oarren-type lathe bed are shown in 
fig. 2.1. The following parameters are taken as the design 
variables in the present work r. 

=» thickness of the main members 

Xj * thickness of the lacing diagonals 

X, = stiffener width on the main and lacing diagonals 
5 

X^ = stiffener depth on the main and lacing diagonals 
Xr- = width of the lathe bed 

5 
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2 .4 .2 Milling Machine Structure 

The structural details of a horizontal knee-tj^^e 
milling machine are shown in Pig.. 2.2. In this case, the followirig 
q’:antities are considered as the design variables t 

= column depth at the bottom 

Xg “ thickness of the over-arm 

X, = width of the machine 
t ■ 

X^ = thickness of the column and table 
X^ = column depth at the top 

Xg = square cross-sectional dimensions of the ribs 
on the over-arm and its joint with the column. 


2.5 OPIBIEJITION HIOBLEM 


Tho complete optimization problem can now bu stated as 


follows . 

i) For lathe bed 
Minimize 


w 

f(S)- 


V. 

D 




% 

1=^ 




subject to 

d^«^ - d^ > 0 

“1 “ ° 


(2.5) 


(2.4) 

(2.5) 


and 


Wg - (co^ + lOO) 1 0 


( 2 . 6 ) 
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< X. < 


a - a 




■i) For milling machine 


Minimize 


f(x) = V. p . + T p 


( 2 . 8 ) 


subject to 


- d„ > 0 

c c — 


( 2 . 9 ) 


0 .'^ - _> 0 


“ (u) + a oo) ^ 0 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


< X. < x(“) 

a - a - a 


a = 1 ,2, . . .,n 


( 2 . 13 ) 


¥fhere 


f(^) = v/eight of the structure (objective function) 

"blri 

V. = volume of the j finite element (the area and the 

a 

thickness of the plate are variables in, case of 
plate elements, and the (square) cross-sectional area 
and the length are variables in case of frame elements) 

p . = density of j^^ element 

J , 

X. = j'*'^ design variable (superscript 'L’ and 'U' indicate 

D 

the lower and the upper limits on X ) 

J 


= number of plate elements 
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Wp = number of ftame elements in which the design 
variables are involved, 

d. = angle of twist of one end of lathe bed when the 



other end is fixed 


= upper bound on the angle of ti-vist of one end of 


lathe bed when the other end is fixed 


d^ = maximum deflection of the cutter centre- (milling 
’.machine) in any direction 



G 

MIN 




tu. 


1 


= upper bound on the maxiraum deflection of cutter 
centre (milling machine) in any direction 
= minimuin negative in phase cross rooeptance of the 
cutter centre relative to the table (milling machine) 
= upper bound on the minimum negative in phase cross 
recoptance of the cutter centre relative to the 
tabl.e (milling machine) 

= lower bound on tho first natural frequency of 

vibration 

th 

= j natural frequency of vibration 


Here the quantities d^^, d^ , are implicitly 

dependent on the design variables and therefore the constraint 
equations (2.4) to (2.6) and (2.9) to (2.12) are called the 
behavior constraints. The method of evaluating tho behavior 
quantities is discussed in chapters 4 and 5. Equations (2.7) 
and (2.13) represent the geometrical or side constraints, which 
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itiipoao limitations on the size of the design variahlos. The 
geometrical constraint equations are linear where as the helmvior 
constraints and the ohjective function are nonlinear and hence 
the prohlsm formulated is a nonlinear mathematical programming 
problem. 



memDers F 11 Main members 





FIG-2.1 WARREN-TYPE LATHE BED DETAILS 
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CmPTER 3 


IDEILIZATION OP MACHHffi TOOL STRUCTOEES 

SrbGnsivo efforts bare Ijeon expended over the years for 
tho development of techniques for the static and dynamic analysis 
of complex structures. With the advent of modem high-speed 
diprital computers, the find to clement displacement method has been 
flovolopod as one of the most convenient tools for analyzing complex 
stmactun.'S. There are several v/ays of modelling the structure by 
tho finite clement method. Therefore, in this work, a comparative 
study of the static and the eigen value solutions of Warren beam 
and milling machine structures is made hy using two different 
triangular plate bending elements. For both the elements, the in- 
plane displacement components u and v_ are assumed to vary linearly 
as • 

(3.1) 

(5.2) 

whom u. and v indicate the components of displacement along tho 
local X and y axes as shown in Fig. 5-1 • ^or the first clement, 
a 9~tGrm transverse displacement model, valid for the vfhole 
triangular element, is used.^ ’’ ^ The displacement model is 

given hy 


u(x,y) :: + Ug X + a^ y 

v(x,y) - + a^ X + ag y 


i 
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w(x,y) - 4 a^x + a^y + a^x + + ag xy + x^ 

2 2 ^ 

+ aj^(x y + X y ) + a^y^ (5.5) 

where w(x,y) represents the component of displacement along a 

direction transverse to the plane of the plate as shown in Fig.3.1. 

In the case of the second, element , the triangular element is divided 

into three parts by joining its centroid to the three vertices as 

shown in Fig. 3-2. A oix-torm transverse displacement model is used 

42 

for o.?,ch of the sub-triangles as s 

''„(x,y) « + ^2^ X + y + i * sy, 

m - 1,2,3 (3.4) 

whore '’f^jn(x,y) is the component of displacement of the sub- 
triangle alon'' the local z-direction . 

In three dimensional structural analysis, both the in-planc 
and the bending stiffnesses have to be considered. The in- plane 
and tho bending stiffnesses are uncoupled for small displaceiaents. 
Those stiffnesses are combined to form the elemental stiffness 
matrix in local coordinate system as shown in Appendix A. Owing to 
the difficulty encountered during transfomation, 62 is not coming 
into tho analysis As shown in Appendix A, and its conjugate 

force M„ are included in tho -inalysis by inserting appropriate number 

Zt 

of zeroes into tho stiffness matrix. The following section gives the 
numerical results given by tho two finite element idealizations. 
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in roforonco 5 • Even in this structure it can he observed that 
the results obtained using the displacement model given in 
equation 5*5 closer to the values reported in reference 5» 
The computer time required is also considerably smaller for this 
moiel as given in tabic 3»5« Hence this model was selected for 
this work. An important thing to be noted in the present study 
is that a lessor number of finite elements yielded results that 
are corapnrablo with those obtained by using a larger number of 
elements in reference 5 . 
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MAETSIS OF CUTIMG POEGFS MU CHaTTEE STaBILITT 

Kie static aad dynami*. rigidities of machine tool structures 
nave an important influence upon both surface finish and metal 
removal rates. In this chapter, the force analysis in horizontal 
milling machines and lathe machines is made in order to evaluate 
static rigidity constraint. A method of including regenerative 
chatter stability constraint in the design of horizontal milling 
machines is also developed. 

4 . 1 HCEIZC:: -AL i'lILIING MACHINE 

4.1.1 Chip Formation in Plain Milling Process 

In plain milling operation, the milling cutter removes chips 

of var.ving thickness intermittently. 'The geometry of xhe plain 

-45-50 

milling process has been studied in a number of works . A 

48,49 

detailed analysis has been reported by Martelloti ’ who clearly 
showed that the path of a milling cutter tooth is trochoidal and 
therefore the milling chip cross sections are confined within the 
arcs of two trochoids. As shown in Fig. 4.1, tiae cutter rotates 

against the direction of the feed in ’up’ or, conventional milling, 
whereas, the cutter rotates along the direction of the feed in ’down 
milling'. Martelloti has developed the following equation for the 
maximum uirdeformed chip thickness (t^) : 
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f 


.J, - 

m 




V.d 


cos 6 


(4.1) 


where, 

D = cutter diameter, mm 
d = depth of cut, mm 

V = feed rate, mm/minute 

V = cutting Telocity, mm/minute 
e = helix angle, degrees 

z = number of teeth in cutter 
s. = feed, mm/tooth 

u 


In equation 4.1, the upper sign refers to 'up' milling and 

the lower one refers to 'down’ milling. If V » t, the trocho- 

51 

idal arcs can be approximately taken as circular arcs . Ihe 

assumption of circular path for the tooth implies that the feed of 

the work piece is inteimittent, so that the woik piece is stationaiy 

during the engagement of a tooth and then mores suddenly by a 

57 

distance of s, before the next tooth starts cutting . This 
assumption is aecessaiy to arrive at simple approximate relations# 

Further, if (^) « 1, equation 4.1 can be written as; 

t s , cos 0 (4.2) 

m ^ t D 

From Fig. 4.2, 


2 d 

cos 4 ) = 1 - — 


(4.5) 
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Hence, 

sinv;^2/^ (4,4) 

Therefore equation 4.2 can be written as 

®t' ijj . cos 6 (4.5) 

Refering to Fig. 4.3, the uncut chip thickness at any instant is 
given by: 


t . = s . sin ifi 

J- 0 1 

The average uncut area of cross-section traversed by one tooth for 
a strai^t teeth cutter is given by: 


A = B s sin (^) 
avg t '2 


(4.6) 


where B is the width of engagement 

Wien more than one tooth is in contact as shown Big. 4#4j the 
maximum uncut chip area is given by: 


z, 

1 


max 


== B s , y sin lb. 


(4.7) 


i=1 


where is the number of teeth in contact at anj' instant. 


The average cross-section of the chip for helical and 
straight teeth cutters when several teeth are in engagement 
can be calculated as follows. 
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B4 


Bati material removed 


s. z = irD A 
t avg 


in one revolution of the cutter = 


Hence, 


Bd s z 

A = 

avg ttD 

where A is the average cross-section of the chip due to 
avg 

number of teeth in engagement. 


(4.8) 



4.1.2 Ibrce System in a llain Milling Process 


Figure 4.5 illustrates the milling action of a strai^t teeth 

cutter simulated by a rotating single point cutting tool. The 
54 

Merchant's circle diagram for cutting forces is also indicated 
in this figure. Figure 4.5a gives, for 'up' milling, 


F 


■R 


and R 


where, 


R cos (?. -y), 


R sin (X-y) 


F 

s 


cos (b+X-y) 


Fj = tangential force on the cutter 

F = radial force on the cutter 
R 

R = resultant force on the cutter 

F = shuar force on the shear plane 
s 

X = friction angle 
y = normal rake angle 
^ = shear angle 


(4.9a) 


(4.9b) 

(4.9c) 
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Thus 

F . cos (X-y) 

p = 

cos (qj+A-y) 

F . sin (X-y) 

and F^ = ^ 

■«os (4 + X-y) 


(4.10) 


( 4 , 11 ) 


The shear force F can be expressed as 
s 


i-^ = B sin (t.) 


( 4 . 12 ) 


The substitutijdg of equatioa-4* 12 into equations 4.10 and 4.11 
gives 


F 


T 


B T , s. . sin 'll. cos (X-y) 
S t 1 

sin 4 • cos (i+X-y) 


( 4 . 13 ) 


and 


B T .s,. sin 4 '-. sin (A-y) 

3 ti 1 

sin . cos ( 4 +X-y) 


(4.14.) 


The horizontal 'and -vertical force components acting on the cutter 
can be expressed as: 


Fjj = Fj cos 4 . + Fjj sin 


B . X .s. sin 4 . 
s ,t 1 

sin 4 . cos (4+X-y) 


[cos i'..cos (X-y) + sin 4. .sin (X-y)] 
1 ’ ^ 


B. 1 '. s,. sin 4. 
s t 1 

sin 4 .cos (d+X-y ) 


COS 


(x-y- 4^) 


(4.15) 
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aad 

sin - Fj sin 
B. T . s . sin il-'. 

® F i_ j- cog sin (A-y) - sin cos (X-y)] 

sin 4> cos (dj+A-y) 


B."!; .s . sin 

S u X 

sin <{) .cos ((^+X-y) 


sin (X-Y- i'^) 


(4.16) 


ViTnen more than one tooth is in contact, equations 4.15? 4.15 and 
4.16 can be generalized as 


B.T .s. cos (A-y) 
s t 

sin <j) .cos (4 +4- y) 
B.t . s, 

s t 

sin 4 .cos (4 +4- y) 


B.T . s. 

St 

Jl _ — 

sin 4 . cos (4+4 -y) 


(4.17) 


'l sin cos (A— Y— 4 j^) (4.18a) 

1 

^ sin 4 '^ sin (A— y- 4 j^) ( 4 . 18 b) 

1 ^ 


sin 


where z is ihe number of teeth in contact at any instant 
i 

and T is the dynamic shear stress. Similarly, for down 
s 

milling, it follows that (from Fig. 4.5b), 

j, ^ — ^ sin 4 ^ cos (A— y+ 4'^) (4. i9a) 

^ sin 4 -cos (4 +A-yO 


and 
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P = ^ 1 sin ip. sin (X-y + ik ) (4.19b) 

sin 4 . cos (b+X-y) i 

It is evident from equation 4.18 and also from experiment 1 
results obtained by various authors^'^’^^’'^^ that the vertical 

component of the force is very small in 'up' milling and its intensity 
tends to decrease as the depth of cut is increased. A reversal of 
the direction of the vertical component of the force may occur when 
the depth of cut is large. In down milling, there will be no 
reversal of the vertical of the vertical component of force as the • 
depth of cut is increased. 

In 'up ' milling for practical values of X , y and d the 
vertical force component will usually fluctuate between -0,21^, 
and +0.2 The positive sign corresponds to the vertically 

upward direction and the negative sign corresponds to the vertically 
downward direction. In both the methods of milling, horizontal 
force components of considerable magnitude will occur as reported 
in references 47, 48, 49 and 55- This behavior is also evident 
from equations 4.18 and 4.19. In .'up' milling the usual range 
of Pg is given by 0.8 to 1.0 P^. 

The variation of the chip cross-sectional area and, hence, 
the cutting forces in plain milling operation with a helical cutter 
are dependent on the values of B,6, d and where 4^ is the 

angle by which the front and the back cutting edges of a helical 



45 


52 

tooth differ. Opitz has shown tnat the individual tooth pulse 

can be approximated by a sawf-tooth shaped curve when ^ j 

53 

and by tra^pezoidal curve when 6^> ip . Singlet, al. have 
experimentally observed such a periodicity in the cutting forces, 
’vhen , where m is an integer and is the pitch 

angle at the same cutting edge, the resultant tooth pulse would 
be a smooth curve. If m is an integer, the tooth pulse would 
be given by a constant superimposed by a pulse variation. 

55 

Boston et. al. have experimentally measured the cutting 
forces in plain milling operation with helical cutters. With 25° 
helix angle, the axial force (B^. ) was found to be approximately 
of the tangential force for various tahe angles. 

56 

Vulf has shown that the average during plain milling 
can be estimated as; 


■T''avg 


(t f 

^ avg" 


A 

avg 


where, 



mean chip thickness in plain milling, i.ni 

s^. sin 77 . cos 6 
t -2 


(4.20) 


and c, k are the constants depending on the material of the woifc 


piece and the environment. substituting the expressions of 

b 

t , sin— and A , the eouation 4.20 can be written as 
avg ’2 avg ’ 



B 


1 


. z 


(4.21) 





(cos 0 ) 


( 3 ^) 


1-k 


(-) 


In the present woifc, the static forces acting on horizontal 
Eiilling iiiachines are calculated using the following parameters: 


D = 10Qnii:i 
3 = 90 nun 

z = '^2 teeth 

s^ = 0.1 nmi/tooth 

ij; = 30° 
e = 25° 

c = 140 (for mild steel) 

k = 0.28 (for mild steel). 


The orientation of the cutting forces in 'up ' milling are generally 
more unfavourable from the point of view of the static ridity of 
cutter centre and also from chatter stability point of view. Hence, 
in the present work, only 'up'' milling is considered in the desigii 
of horizontal milling u^ichines with the following values of the 
cutting forces: 


H 


X ‘avg 


''t ' larg' '"S- 

0.2 168 kg. 

= / 

H V k 


R 


= 875 kg. 



47 


4.2 lOECE SlSmi IN lATEE DURING TURNDIG 

!I!he cutting forces in the general case of turning operation 

51 

are shown in Dig. 4.6. (Die resultant force E is given by 



where 


P ,= feed force on the cutting tool in the direction 

JCi 

of the tool travel 

P^= thrust force on the cutting tool in the direction 
perpendicular to the generated surface 
P^= cutting force on the cutting tool acting in the 
direction of the cutting velocity vector 


In view of the fact that the twist of the lathe bed fortis a large 

part of the total deformation, it will have a considerable influence 

58 

upon the wording accuracy . Hence, torsional rigidity is included 
as one of the behavior constraints in the design of lathe beds. In' 
lathe machine, the work piece and its support and the cutting tool 
and its support are much more flexible compared to the lathe bed. 
Hence the contribution of lathe bed for chatter in turning operation 
is not sigiiificant compared to the other machine members indicated. 
This is the reason why chatter stability is not included as a 
constraint in the design of lathe beds. ' The maximum torque (t_. ) 

on the work piece is calculated as: 



4i 


2 frU . r- = H.I. 
min max max 

where U . is the minimum spindle speed, and H.P. is the 
min ^ ^ max 

maximum horse power of the lathe. 

The twisting moment acting on the lathe bed depends on the 

work piece diameter and the orientation of the cutting forces. 

In actual machine, it is very difficult to estimate the torque to 

which the bed is subjected. Hence, T is taken as the twisting 

max 

moment acting at one end of the lathe bed (considering the other 
end is fixed) to compute the torsional deflection. 

4.3 GiSK'iJRAl PMiaRjilS OH CHilTER IH MAOHIIE ODOIS 

Severe vibrations in metal cutting machine tools are often 
encountered giving rise to undulations on the machined surface and 
excessive variations in the cuxting forces. Such energetic vibra- 
tions belong to the class of self-excited vibrations, commonly 
known us 'chatter’ in machine tools, and the source of self excited 
energy is derived from the cutting process itself. As presented 
in Pig. 4.7, the self excited vibrations are due to the interaction 
of cuitiug forces and the machine tool structural dynamics. Seveia,! 

researchers have proposed various models for the self excited 

59 60 

vibrations in machine tools. Kasirin and Arnold have utilized 
•vs.'locity principle' for explaining chatter in machine tools. !Ehe 
velocity pxlnciple assumes that the characteristic of the cutting 
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force is that it has a negative slope with respect to tne velocity 

of cutting. Ine regenerative chatter arises because of the 

variation of the cutting forces occuring due to the variation of 

the uncut chip thickness j and it is explained with respect to Mg. 4.8. 

Mr the sake of brevity, only a single degree of freedom ^stem is 

illustrated in this figure. The cutting tool cuts a surface which 

it hj,s alreody cut during the previous revolution in case of turning, 

during the previous stroke in case of plaining and by the preceding 

15 

tooth in the case of milling . Therefore, from Mg. 4.8, if there 
was a vibration betw'wen the tool and the work piece during the 
previous ith cut resulting in an undulated machined surface, the 
chip during the next (i+1 ) cut is removed from the same undulated 
surface. Hence, the amplitude of variation in chip thickness, AY, 
is given by 

AY = Y. , - Y. (4.22) 

1+1 1 

where 

Y^ = vibration amplitude betv/een tool and woifc piece 
noimal to the cut surface in ith cut 
Y^_i_^= vibration amplitude between tool aM work piece 
normal to the cut surface in (i+l)th cut 
If ¥ is the phase difference between Y^ and j then AY 

can be written as 

ay = qY^e“^’- (4.23) 


v^fhuro <1 is a constant. 
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Ihe change in. the cuxting force (hP) due to the variation 
in the uncut chip thickness (ay) is given by 

Ap = - B.r. AY (4.24) 

where r is a coupling constant depending on the cutting conditions, 
and B is the width of chip. OZhic change in the force component 
will again result in a vibration giving an undulated surface on the 
work piece and the regeneration of the undulation proceeds in subse- 
quent cuts. 

61 

Tlusty and Polacek proposed a model according to which 
different modes of vibration of the stmcture of a machine tool 
will be coupled with the ’cutting process dynamics' through cross- 
receptance and not by direct reeeptance. ' According to this model, 
chatter will occur when the energy delivered to the vibratory system 
by the cutting forces exceeds the energy dissipated by the vibratory 
system. Ihis model is known as the 'mode coupling theory*. 

In actual practice, all these three phenomena appear simul- 
taneously and out of the three theories it is generally considered 
that the generated energy as given by the velocity component principle 
is low compared to those given by the regenerative chatter ai'jd mode 
coupling principles'^. Tlusty and Polacek have taken the simplest 
dynamic cutting-force relationship by assuming a direct proportionality 
between the force and the undefoimed chip thickness. In the analysis 
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H. , the fore, relationship mn 

made by Tobias, et. al. j 

1 . the rate apart from the 

related to the cutting speed —-i ™ 

undeformed chip thichness. 

TO derive the limits of stability for the machine 

tho cutting s.ste., it is neosssat. ^ assume a iorca relationship 

this to the machine dynamics, 
for the cutting process and relate thi 

rtf he solved directly, but 
The governing equations of motion cann 

V. V\^T oqg nTr>i, ng ijb.© condi*ti-^^ 

the analysis can be simplified very mu 
at the limit of stability. 

4.4 mathematicai expbessiou for the limit of stab 

i. ■ thsaorv of stability analysis 
In the present work, the basic 

■t, n 1.15,61 T;ised in applyi-i^S etiuation 

as given by Tiusty and Polaceh 

-HiT-ptr The simplifyf^S 
4,24 for the design of milling machine s ruco 

assumptions made are: 

a) The titration ^stem of the mschlne is imear 

b) The variable component of the cutting f 

direction (y) remains constant 

+ n-p the cutting force depends 
c) The variable component of t 

• +h^ d irection of the normal to 

only on vibrations m th 

the cut surfQ-ce 

4 ) The value of me variable component of the cutting 

force (41) varies proportxonat ly LIBRARY 

vdth variable oisplacatent, ^ _ A 
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e) Ihe frequency of vibration and the mutual phase shift 
of undulations are unaltered in subsequent overlapping 
cuts. 

Equation 4.24 can be rewritten, in view of equation 4.22, 
as 


AP = 


Br (Y - Y . ) 
^ i+i i^ 


(4,25) 


At the limit of stability, the variable component of cutting 
force Ap will result in cross -re ceptance whose .-amplitude 

is same as Y^. Hence 


Y . = AP. 4 (m ) 


(4.26) 


and 


Y. 

1 



where $(u)) 
sy stem 


= 1 ■ (4.27) 

is the absolute cross-receptance of the vibratory 


(oi) = 5(a)) + jH(o)) 

where, G(m) and H(a)) are the real and imaginary components of 
4(a)) respectively. 

Here 4 is taken as a particular case of cross-receptance in which 
the variable cutting force ( Ap) is along the direction of the 
cutting force and the resulting relative vibration (Ay) between the 
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tool and the vrork piece is along the aomaal to the cut surface. 
Combining the equations 4.25, 4.26 and 4-27, the following 
equation is obtained: 


l/(Br) + G + jH 


G + jH 


= 1 


(4.28) 


Since the imaginniy parts of tne numerator and the denoih-nator of 
equation 4.28 are equal, it follows that, 


l/(Br) + G 


Equation 4.29 is satisfied only if, 


Br 


+ G = -G 


(4.29) 


or 


JL 

2Br 


G 


(4.30) 


For certain specified cutting conditions, tbe value of the 
coupling coefficient (r) is fixed^ Therefore, the liroiting 
v.jlue of chip width "fehe limit of stability, for these 

specified cutting conditions, is governed by the minimum negative 
real cross-receptance ( - G^.^^^ ) of the structural system. Tfaere-^ 
fore, the lowest limit of stability is deteimined fraa the 
following equation 




- G. 


LEi 


(4. 31) 
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The meaning of equation 4.31 can also be explained with reference 
to li-g. 4.9. 

4.5 DETEfflIMTION OE GOUIEIEG CONSTMT IS PLAIN fflILIING 


Equation 4.31 is used, for the design of horizontal milling 
machines for regenerative chatter stability, in the present work. 
The value of coupling coefficient (r) depends on the cutting 
conditions. If the critical or the maximum value of coupling- 
coefficient estimated from all the commonly used 

cutting conditions, structure can be designed for chatter stability 
by putting a limit on the negative inphase cross receptance 


(a 


LIM 


^^1114 ^critical 


). Therefore it becomes essential -to 


analyze the various cut-feing conditions inorder to estimate r 

critical 

so that -the value of • can be deteimined, 

JjIM. 


4.5.1 Coupling Coeffiqient for Strai^t Teeth Cuttt=r from Basic 
Theory of Metal Cutting. 


In Pig, 4.10, AY and AP are shown for a strai^t teeth 

cutter in horizontal milling machine, Refering to Pig. 4.4, 

equations 4.18a and 4.18b for ’up' milling can be v/ritten as 

r, z. 

B Tg 1 

Eg I (s sin 4 ). cos (l-Y- -4. ) 

sin 4 cos(4 -tA-y) 1 


(4.32a) 


and 
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B T 


p ' = 


sin (J) cos(({)+^- y) 1 


y (s^ sin Ti,. ). sin (x-y-4) 

^ U ^ 1 i 


(4.32b) 

where s, sin indicates the uncut chip thickness at 1th 

tooth. Refering to Rig. 4.10a, for a relative displacement 

(^y) of the cutter centre with respect to the v/ork piece, the 

change in uncut chip thickness at i^^ tooth is given by aY.cos 

th 

where ^ . is thu angle between i . tooth and the nonaal 

J 

direction to the cut surface .tS shown in Eig# 4.11« 


* 


The corresponding new horizontal force and the vertical force 


become 


* 


B. T 


sin 4 cos{<(’+X^y) 


z. 

1 


I (s^.sin cos cos (X ) 


and 


s 

7 


(4.33a) 


z. 

1 


y (s, .sin t|; -ay. cos Ip . ).sin (X-Y-^p. ) 

sinbcos(v+X-Y) ^ 1 ^ 


(4.33b) 


The changes in the forces ure given by 


AP = 1 - P = 


B. T 


1 


y Alices * cos(X-Y*^ . ) 

sin3cosOP+X-Y) \ 


(4.34a) 


* 

A P = P “^P ^ 


B.T 


2. 

X 


sin <!> cos(b + X-Y) 1 


y AY.cos ip ..sin(X-^-ijj ) 




(4.34b) 
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Using equation 4.24, the coupling constants can he expressed as: 


AR 


H 


z. 

1 


3. AY 


sin 4 cos('H^y) 1 


I cos cos 


(4.35a) 


and 


AR, 


T B.AY 


^ — : I cos R. sin (X-r-'j'. ) 

sin d cos(4 +X - t) 1 ^ 


(4.35b) 


where r^. and r^ are the coupling constants in horizontal and 
vertical directions respectively. 

Therefore 


/ 2 2 

rg + 


(4.36) 


Siiriilarly for down milling, 




sin 4 cos((J<+X-y) 1 


1 cos tR . cos (>-y+i|>^. ) (4.37a) 

t/ 


and 


Ty - 


z. 

1 


2 J 

sia4 cos(d+X-Y) 


cos sin (X-r+il)^) (4.37b) 


In the present work, the relationship due to Lee and 
64 

Shaffer is used for coiaputing the shear angle; 


‘■ = J + Y-X 
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Por ^ = 18° and y = 10°, this equation gives ^ = 37°. 

. 65 

The dynamic shear stress, t , is computed from the Abuladae's 

s 

shear stress relationship for ductile materials: 

— 6e 

= 0.74 a (6) (4-38) 

o U. 

whore is the ultimate tensile strength of the material and 

e is the percentage elongation. 

2 

Por = 40 kg/mm , e = .2 the value of obtained is 

36 kg/mm^. 

The values of coupling coefficients (r) obtained by using 
equations 4.35 to 4.37 for various values of ^ are given in 
Table 5.1. The cutter teeth are taken to be symmetrically disposed 
to the mean direction Y(noimal to the cut surface). 

4.5.2 Coupling Coefficient from Vulf’s Formula 

Eefering equation 4.20, 

p I S A 

■‘^T'avg “ / vk ' avg 

^’^avg'^ 

Since Jg = , IV - 0-2 Ijlavg > 

R =/ljj + ~ ^lUvg 

E 1 
T'avg 
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Hbnce^ 


R = — r . A 

(t f 

avg^ 


(4.39) 


If the displacement of the cutter along the noimal direction 
to the cut surface is ^‘*'a_yg corresponding diange 

in the cutting force is AE^. 


As Ar 0, one obtains from equation 4.39j 


(4.40) 


7 ? =-Br=c. A . (^) (t )' 

^^ayg ^\yg 


( 4 , 41 ) 


Substituting the value of A.^yg f^xim equation 4.-8, equation 
4.41 can be written as 


r = c (|) . (-|) z.k 


(4.42) 


since t = s, . sin — • cos 6 ,■ 

LlTg t 2 ’ 


c. z.k 


(-) ^ 


7T(cOSef'"' (S^r 


(4.43) 


By substituting the values of c,z,k, 6 and s^ given ii 


section 4.1.2, in equation 4.43 »• 


r = (331) (d/D)* 


(4.44) 
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The value of for various values of (■^) are given in Taole 
5-2 along vdth the corresponding values of . 

It can be observed, from the tables 5.1 and 5. 2, that the 

values of r obtained by using equation 4.44 agree veiy well with 

those obtained by using the basic theoiy of metal cutting. In the 

present woik the critical value of (r ) is taken as 165 kg/mm *. 

LIM 

The width of cut (b) is taken as 9 cm. 

In the chatter stability constraint, the upper lii'ait on 
is taken as 

Ic^^h - . 1 

' mu' ~ 2 B r^^, 

LIM 


= .f)00003 cm. /kg. 
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TABIIB 4.1 

YAIUES OF GOUPUKfO COFSTaHTS EROM BaSIG THBOHT OF 
WEHAL GUTTING 



■ ■ ‘ -“2 

Goupling constants r (kg/mm ) 


'up' milling 

'dovm' milling 

350 

179 

158 

40* 

159 

157 

45* 

159 

156 


TiBIE4 .2 

V aides of GOUPtrUG GONSTaUTS 

FROM avIPIRIGAl Fomcrial 


d/D 

2 

r (kg/mm ) 

35° 

.09 

139 

40° 

.116 

152 

e 

45 

.146 

165 




ARC^ef^lS A TROCHOID 


(a) UP MILLING {b)DOWN MILLING 


FIG. 41 THE PLAIN MILLING PROCESS WITH STRAIGHT 



FIG. 4-2 GEOMETRY OF THE PLAIN MILLING PROCESS 





FrG.4-3 THE UNCUT AREA OF CHIP-CROSS SECTION 

(CONFINED WITHIN TWO CIRCULAR ARCS) TRAVERSED 

BY ONE TOOTH 





(3) UP MILLING 


CUTTER 




WORK PIECE 



(bj DOWN 

MERCHANfS CIRCLE ni&fiRAM FOR CUTTING FORCES 
^ A ROTATING 


'•:, i'; ‘ 


ornMT r-i ITTIWA TOni 



FIG.A-6 niTTING FORCi^S IN CONVCNTION&L TURNING 






FfG»47 BASIC DIAGRAM OF CHATTER IN MACHINE TOOLS 



y- indicates the drection of the 
normal to the cut surface 


fig4'8 basic diagram representation for the process 

OF SELF EXCITED VBRATIONS IN MACHINE TOOLS* 



€S 


G 



FIG,4-9 REA L CROSS- RECEP TANCE CURVE WITH RESPECT 

TO FREQUENCY 




FfG..4'10 ORIENTATtON gUTiNG FOf^ES AND THE DIRECTtC^ 

OF THE NOR MAL TO THE CUT SURFACE iH 
' MILUNG : ' ' ■ 



FIG.4 ir ORIENTATION OF GUTTING TEETH IN ACTION 


WITH RESPECT TO THE NORMAL TO THE 


CUT SURFACE 
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S'i'irPlC AHD DYIMIC MllYSIS OP MCHn® TOOL STRuCTCJELS 

■5.1 LPPLECTION AM) STHESS illAXYSIB 

Once the load vector P, for ’.vhich the machine tool 
stricture is to "be designed, is determined as explained in the 
previous ch,;i,pter, the displacement vector Y can he computed by 
solving the following equations ! 

[k] Y = P (5.1 • 

vhQTc [K] is the master stiffness matrix of the structure. 

In the present work, the solution of equation 5*1 is obtained 
by the oholc?sky decomposition of the symmetric banded matrix, 
storing only the upper triangular matrix, followed by forward 
and biickward elimination technique. The horizontal milling 
machine structure shown in Pig. 5*1 is analysed for static 
deflections and stresses. The idealisation of the structure 
is sho\m in Pig. 5*2. The dove-tail (vortical guide for the 
knee) thickness on the front face of the column is talcon as 

5.5 cm. The dove-tail thickness on the over-arm is taken as 

2.5 cm. These thicknesses are added to tho thickness of the froni 
face of the column and the hcttoia side of the over-arm respectii 
The arhor support is idealised as two frame elements. The 
rectangular cross-section of these frame elements are taken as 

7 X 10 cm. The arhor diameter is taken as 4 cm. 



The maximum deflection of the cutter centre in any 
direction is .d066cm.In metal cutting imchine tools, the induce 
stresses, corresponding to the permissible deformations, are 
generally far less in value than the allowahlo stresses for 
various materials. For the structure shown in Fig. 5,1, the 

maximum principal stress developed in the structxire is 11.97 Kg 

2 

and the maximum shear stress induct-d is 6.63 Kg/cm . It is ohv 
that these values are far helow the allowable ones. Hence in t 
design of machine tool structures strength was net taken as a 
design constraint. The stresses are computed from the known no 
disulacements hy using the stress-strain and the stra.in-displac 
relations of linear elasticity. 

In the design of horizontal milling machini-s, the 
maximum cutter centre deflection should not exceed the upper he 
d^,^^ and this requirement can he stated in the optimum design 
problem as 

d^'^^-d^lO (5.2 

c c 

EIGfillS VxlhTJE EROELFIi 

In the finite element method, the eigenvalue problem 
solved by choosing the nodal displacements as generalized ooorc 
mtes. Thus, the linear eigenvalue problem given by the folloif 
equation 

IK] Y = AM Y 


(5.; 



Pig. 5 . 3 . Since the eigen vectors are required to be ortho-noima 

<y 

the reduction of the generalized mass matrix [M ] to an identily 
matrix is used as a convergence criterian in sub-<space iteration. 
Por the milling machine structure (Pig. 5.l), the natural £reque 
and the mode shapes are shown in Pig. 5.4. Por the lathe bed (Pi 
the natural frequencies and the mode shapes are shown in Pig. 5.6 

5.3 PYIAMIO EPSPOSSE OP IHE SIMJCODBE 

Ihe dynamic response of the structure has been obtained by 
using modal coordinate system* The damping matrix [C ] is taken 
proportional to a linear combination of the stiffness and mass 
matrices. Since it is not possible to estimate the modal damping 
factors from the blue prints of a given structince, the values ol 
damping factors have to be obtained from experimental results on 
similar structures. In the present work, in the desd@3. of horizc 
milling machines, the damping factors are assumed to have a value 
0.06 for the first few modes. choosing the nodal displacemeat 
as generalized coordinates, the differential equations of motion 
structure can be written in the following foimj : 

^ ^ ?■ 

[m] Y + [c] Y + [f] Y = P sin mt (5*4 

where, P is the vector of force amplitudes acting at the various 
nodes of the structure and m is the forcing frequency* 

py using coordinate transformation 


Y = [U] n 


(5*5 
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in v/hich [U] is the [M] -orthonormal modal matrix, and it 
represents the vector of modal or normal coordinates. 

Si'om equation 5.5, 

i = [U] I (5.6) 

Y = [U] n (5.7) 

Substituting equations 5.5, 5.6 and 5.7 in equation 5-4, and by 

1 

premultiplying throughout by [U ] , 

[Ur [M] [u] n + [U]^ [C] [U] n + [U]^ [K] [U] n 

= [U]^ f.sin a)t (5.8) 

T 

where [U] is the transpose of [U]. Since the modal matrix is 
orthonormal with respect to the mass matrix [i|i , and the damping 
matrix is proportional to a linear combination of [K] and [M] , 
equation 5.8 results in a set of decoupled differential equations ; 

[“•I-.] n + ['"2^'^-] n + n = q sin mt (5.9) 

->■ !I! 

where, Q is equal to [U] j and C represents damping factors. 

If only p natural modes out of JJ are used in [U] , equations 
5.9 represent p uncoupled second order ordinary differential 
equations. Ihe equations 5.9 can also be written as 

+ 2?^ 0)^ + 0)^ nj_ = sin wt (5.10) 

Equation 5.10 can be solved to find the normal displacement n.(t). 
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(These displacements, n^(t), can be used to obtain the dynamic 
displacement amplitude vector, Y with the help of equation 5.5* 

The inphase component G(to) of the absolute cross-re^eptance 
of the horizontal milling machine, shown in lig. 5.1, are 
plotted in Pig. 5.7. In Pig. 5.7, the values of g(o'-) are shown 
in X, Y and Z directions for the cutter centre and the table centre 
separately. Ihe harmonic response loci for Hie same structure 
are shown in Pig. 5.8. 

In the design of horizontal milling machine structures for 

chatter stability, the negative inphase cross-receptances of the 

cutter centre relative to the table have been computed at an 

interval of 10 radians/second beginning from the first natural 

frequency to the fourth natural frequency of vibration and from 

these, the minimua valu^ has been found. (Ehe values of 

fr,,.™ have been evaluated for different inclinations of unit 
MIN 

amplitude of harmonic forcing with respect to the horizontal 
direction (e). Ihe values of for various inclinations are 

given in Table 5.1. Prom these results, it c.^n be observed that 
there is no significant difference in the value of 0- r jjj for varioug 
values of 8 encountered in practical 'up ' milling. 

The chatter stability constraint in the optimum desigp. 
of horizontal milling machine structures can be stated as follow^ f 

4iN ' SilN - ° 


(5.11) 
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where 


a 


(u) 

'MIN 


1 

2 B r 
LIM IIM 


Althou^ the receptances of the lailling machine stiucture 
have been computed taking the modal damping factors as 0.06 f a 
study has been made to find the effect of damping factors on the 
receptances. Bor a 10^ change in the damping tactorsj the response 
loci of the structure are compared in Big. 5.9. Brom the fi^re 
it can be observed that values have changed approximately by 

^ 0 . 


TABLE 5.1 VAHJES OB 0^^ BOR DIBBERENT INCLINATIONS 
OB UNIT AMPLI'niDE OB HARMONIC B0RGIN& 


B (degrees ) 


-20 -10 0 +10 +20 


^MIN 

c i.Ag X 10 


5 


0.6721 0.6769 0.6818 0.6721 


0.6657 




/ / NOTE"— All dimensions 

/ in cm- 

ARBOR Df A = 4-0 

COLOUMN THICKNESS::^ 2.8 
OVER-ARM THICKNESS « 2-8 
TABLE THICKNESS =2.8 
SQUARE RIB SECTtON=^3X3 


FlG.5-1 HORIZONTAL KNEE TYPE MILUNG ,JMAaj^ 







[K] = gonoralisod stiffness matrix 

[M] s generalised mass matrix 

[X] = iiiitial a'j, i^roximate vectors 


i5r 


degrees of freodon 


p = number of required first few node^ , 

X’^nr^R = maximum number of subspace iterations allp\;od 
max 

[Q,] =■ Eigen vectors of the x*i^ojectod operators 


Enter With [K ] , [M] , ], H, p, to 

Set ITER = 0 



contd 






Find tliG mm'ber of rows 'KR’ in which 
tho follovdng criterion is satisfied 

. P . 


I [tiij 

j=1 


< .001 , i = 1 , p 


Find tho oigen solution of the projected 


operators 


[K][Q] = 


Find an improved approximation to the 
eigen vectors 


[X]- [Y][Q] 


Is 'M' < 


ss converged 
Return 


5.3 Rayleigh dlitz Suh Space Iteration Algorithm for 
Determining Eigensolutions 
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FIG.5*4 mode shapes OF HORIZONTAL MILLING MACHINE 



(a) WARREN TYPE LATHE BED 


Stiffeners idealized as frame elements 



5 parts at equal intervals 

(b) FINITE ELEMENT GRID ON BOTTOM SURFACE 
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(c) LACING DIAGONAL (^) sipE PLATE 


■ Stiffeners idealized as frame elements 
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(e) FINITE ELEMENT GRID ON TOP SURFACE ; 

Dotted lines represent, finite element grid 


5-5 DETAILS OF IDE ALIZATION OF WARREN : TYPE 
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R6-57 INPHASE CROSS RECEPTANCETS OF CUTTER CENTRE AND 

table Calong x,y,z directions) of milling machine. 





FIG-5-8 HARMONIC RESPONSE LOCI OF CUTTER CENTRE AND 
TABLE (X.TZ DIRECTIONS) OF MiLUN 6 MACHINE- 
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CHAPTER 6 


OPTIMIZATION METHOD 


In chapter 2, the problem of optimum design of machine tool 
structures has beeii formulated as a nonlinear mathematical progra- 
mming problem. !Qae methods of computing static deflections, 
natural frequencies and the receptances, on which restrictions are 
to be inclosed, have been developed in chapter 5, 4 and 5, The 
choice of the method of optimization procedure and its description 
is the topic of this ehapter. 

6.1 CHOICE OP TEIE METHOD 

The three general classes of widely used nonlinear progra— 

G7 

mming methods are as follows: a) Gradient projection method of 

68 go 

Rosen subsequently modified by Goldfarb : Thou^ this method 

works well with linear constraints, its efficiency is considerably 

reduced in the case of nonlinear constraints, b) Feasible directic 

70 

method of Zoutendijk : This method would also present difficultis 
similar to the gradient projection method, c) Penalty function 
method : In these methods, the (•onstrained forraulcitioa is trans- 
formed into a sequence of unconstrained optimizations which can be 
solved without much difficulty. These methods are quite reliable ■ 
and the Sbqubnti.3iL nature of Ihese methods allows a gradual appronc 
to cri-Gicality of constraints. These methods allow coarse 
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antiroxiimtions to tfe lisod during early stages of the optimization. 

procedure, and close approximations during the final stages, fhe 

coinput£.tional timoj however, is more in those methods hocausc the 

ininimization prohlem has to he solved a number of times. To 

reduce the computer time, the penalty function methods alio?/ the 

use of approximate analysis without effecting the acciiracy and 

stability of the procedure . In the present work, the penalty 

71 

function method of Piacco and McCormick is used as it has "been 
found to be quite reliable . 


n .2 FlinnO- McCORMICK ETTIRIOR PERAITY FJNCTIOH METHOD 


In this method, the objective function is augumentod 
with a penalty term consisting of the constraints as shown belo'w? 


HI 

^ = t( ^ I ~~r7 


where 1)(X ,r]^) is called the penalty function, and r^ is called 
penalty par;?.meter in k^^ minimization step. The minimization of 


6-function is performed for a decreasing sequence of r^ so that 


r ^ 

k+1 


( 6 . 2 ) 


The equation 6.1 requires a feasible starting point and, in the 
present work, it is found by trial and error. Since each of the 
designs generated during the optimization process lies inside the 
acceptible design space, the method is classified as interior 
penalty function formulation. 
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6 , 5 DiVIBOH - FLETCHER - POY/ELL FABIABLD ’!BTRIC IRTCOFSIRAIHED 
fAIFIfllZATIOF METHOD 

In the penalty function formlations, since a sequence 
of unconstrained minimizations has to he performed, the selection 
of an efficient method of unconstrained minimization becomes very 
important . All the methods of unconstrained minimization find a 
sequence of improved approximations to the optiniam according to 
the iteration ' 



X. + S. 


1 1 


(6.3) 


where , 

= the design vector corresponding to the minimum 
of 4' -function along the current direction 
X^ = the starting design vector 

T* = the minimizing step length in the direction . 

There are several methods available for finding the search direction 
in equation 6.3. 

In the present work, the Davidon-Fletchar-Powell Variable 
72 

Metric method is used. This method can bo considered as a quasi- 
Hewton algorithm, and is a very powerful general prooc*dure for 
finding a local unconstrained minimum of a function of many variable 
In this method, the i"*^^ search vector in equation 6,3 is 
computed according to the following relation 

\ = -IH.] V<j> (X.) 


(6.4) 
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where ^<!) (X^^) denotes the gradient of the ^’-f'unction at 
and [ F^] is a positive definite- symiifctric matrix. The matrix 
[f^] is updated according to the following procedure 


t , = [H.] + [A. ] + [ B.] 


(6.5) 


where , 


[A.] = 


[B,3 


T* S. S. 
1 1 

S. Y. 

1 1 


T ^ 

[H.] 7. V. [H.] 

- 


( 6 . 6 ) 


(6.7) 


Vi [ly V. 


and 


V. 

1 


= (\+-,) - ^K\) 


( 6 . 8 ) 


The updating of [H^] presorves the symmetric positive definiteness 

of [H. .] which ensaros the stability of the procedure. The 
1+1 

positive definiteness of is influenced only hy the accuracy 


Tv w 

with which ^ is determined. To start with [H^ ] is taken as the 
identity matrix. 


6,4 ONE DB/ENSIONAJi ?ilINDi!IZATION METHOD 

In the present work, the one dimensional minimzation is 
accomplished hy cuhid interpolation method . In this method, the 
algorithm used to compute 't is reapplied until 
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becomes moderately small. If the cosino of the angle between the 
vectors S and V(f) at the minimizing step laagth t is suffi- 
ciently small quantity (e) , i.c. 

T 

S. V({) 

cose = -r — ; < e , (6,9) 

!sl .!V(>1 

the one dimensional minimization is considered to have converged 
to T , In the present woiic, the value of e is taken as O.O5. 
However to reduce the computer time, the number of cubic inter- 
polations is limited to three , 

6,5 IDIITIOHII. COHSIDEHiVTIOlTS ilH) CCiryEHaS'TCE CRI'IERLl 

6.5.1 Initial and Subsequent Values of r^ 

The initial value of r is taken such that^^ 

1.5 (6.10) 

and the subsequent values of rj^ are found by using the following 
relation 

= 0.1 (6.11) 

^k 

6.5.2 Tie starting the M Matrix 

If a function that is to bo minimized is highly distorted, 
it might occur that is positive after a fc’.ir iterations, 

indicating that ^ is not a direction of descent. In such cases 
the remedy is to set [H ^3 -s [H^] and start all over again. 
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In the present work, however, such a situation never occurred 
in all the examples studied. 


6.5*5 Termination of Ilinimization for Each r, 


k 


The minimization of (j)- function for each value of 

is terminated when the predicted porccxitage difforonce between 

the current and the optimal <!> -values is smaller than a small 
74 

quantity , e ; 


T 




24 ^ 


"i- ^^ i 

24>i 


< e 


( 6 . 12 ) 


The value of e used in the present work varied from 0.05 for r^ 
to 0.005 fo^ subsequent r^ . 

6.5.4 Relative and Global Minima 

In the design of horizontal milling machine strictures, 
tY70 completely different starting design vectors are used for the 
sequence of minimizations to see whether any relative minima exist 
in the design space. Two different starting design vectors led to 
the same optimum design. Although, merely on the basis of two trial 
starting designs, it is hard to say that the miniiiuim oTctain'-d is 
the absolute minimum over the design space, finding the similar 
least weight design by starting from two different initial designs 
is at least a pointer in that direction. 
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6.6 COMHITITIOI OF TBE GRADUCM? OF 6 - FUFCTIOI (Vd) 


The va.riahlc metric method, used in the present vroifc, 
requires the gradient of the ^-function for computing every new 
search direction and the one dimensional minimization method 
might require the evaluation of S. "d more than once in every 
search direction. The gradient of the 6 -function is obtained 
hy differentiating the (' -function with respect to the design 
variables as, 


Vd = 


/3X^ 


( 6 . 13 ) 




In the present work, the forward finite difference scheme 
is used for computing the gradient V6 . This scheme requires (n+l) 
evaluations of the d -function where n is the number of design 
variables to compute . Apart from this a few more d ■i.'function 
evaluations are required for each one diiiienoional minimization. 
Since the evaluation of one d -function itself requires a lot of 
computer time, any alternative procedure, which computes the 
gradient faster, would substantially improve the efficiency of the 
optimization program. 


In computing the d -function at any design' vector, a major 
portion of computer time is spent in the evaluation of the stiffness 
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and mass mtricos and the response quantities i.o., static 
deflections, natural frequencies of vibration and the dynamic 
displacements of the structure. In the present work, an alternate 
procoduro 5 which is more efficient than the finite difference 
scheme, is developed for iiie evaluation of the gradient of . 

This method is based on the reanalysis technique of modified 
structures. The reanalysis procedure deals with the best way of 
obtaining the rosponse quantities at a modified design vector 
given the stiffness and mass matrices and the response quantities 
at a Certain design vector. 

75 

6.6.1 lorivativos of Static Displacements 

Consider again the matrix equations 


[k] f = f (6.14) 

where [K] and P depend on the design vector t. By taking the 
derivatives of equation 6.14 with respect to the design variable 

% ’ 

-> 9[K] . SP 

(6.15) 


r--, 3Y 

3X^ " 3Xj_ 


i.e. 


[iH 


3Y 


-V 

3P 


3X, 


ax. 


3 [k3 : 
3X. 


( 6 . 16 ) 


The terms ^ and oh the right side of equation 6.16 
i i 


are computed by finite difference method . Let the design X be 
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“>• -X- 

porturloed to X Toy changing to X^ +Ax. . At this new design 
X , the stiffness matnix [K ] can "bo generated hy the rogular 
procedure. 


Thus 


rail _ [-K*]- [k ] 
8X^ AX^ 


similarly, 




( 6 . 17 ) 


(6.18) 


where P is the load 
Thoreforo equation 6.16 

r-7,3Y p-p 


IK] 


3X.' 

1 


AX, 


-i-* 

vector at X. . 

1 

can ho written as 

fie*] -[iq + 

[ sr. 1 ^ 

1 


( 6 . 19 ) 


On the right hand side of equation 6.19? there will he as many 

vectors as the number of design variables . Therefore 

this problem reduces to solving an equilibrium equation with 'n' 

vectors on the right hand side. Equation 6.19 has boon solved to 
3 Y 

got ^ by the cholcsky docoEiposition of symiactric banded matrioos 
as discussed in the previous chapter. The displacements at X 
computed by using equation 6.19 compared with those obtained 
by the exact raothod in Table 6.1 . It can be observed from this 
table, that the values given by the present approximate method are 
al.most same as those obtained by the exact method. 
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6.6.2 Derivatives of the Natural Frequencies 

Since the first few natural frequencies of vibmtion 
of a gi„neral eigenvalue problem 


[K] f = w^[M] Y 


( 6 . 20 ) 


are only inyolved in the constra,ints, the values of -57 ? 
k 1 required for use in the reanalysis cycle of 

optimization. By differentiating equation 6.20, the exact 
expression for the rates of changes of eigenvalues is given by 


3co 


1 




T 5[k] 


[Y. [ 


I Y 


m. Y. 
1 1 


T 3|l£ij]Y. 

( 6 . 21 ) 


This equation has been used for the eigenvalue analysis at 

21 

modified designs by Fax and Kapoor . The partial derivatives 
of [K] are found by using equation 6.17 and those of [M] using 
a similar equations 


r 9 M . - [m 3 

A X^ 


Thus equation 6 .21 becomes 


( 6 . 22 ) 


9 03. 
1 







2 

03 Y. 

i 1 




t 


2 

03. 

1 


"> T r** 
1 


[m 3 


Y. 3 

1 


( 6 . 25 ) 


The -values of at X calculated hy using equation 6 .23 are 
compared with those obtained hy the exact method in Table 6.1 . 

It can be seen that the values obtained by equation 6.25 are in 
excellent agreement with those obtained by the exact method. 

The values of as predicted by the approximate method 
(equation 6 .23) and by the exact method for various percentage 
changes in the design variable are shovm in Tabic 6.2 and in the 
Pig. 6.1. These results indica,te that the frequency values do 
not differ by more than 2 .'fjo even for a change of 2C;'o in the 
value of the design variable. 


6.6.3 Derivatives of the Eigen Vector 

In the present work; a partial solution to the general 
eigenvalue problem 


[k1 Y = Y 


(6.24) 


is solved for computing the reoeptanoes of Iho structures. 

Since the eigen vectors form a complete sot of vectors, any H-* 
component vector can bo represented as a linear combination of the 


II -eigen vectors as 
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Z ,N 
^^i . y 

cK. ” ^ 


k=1 


a . Y^ 
xjk k 


(6.25) 


where, 
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a. . „ = 
10 I 


M - V -S-’r,, 

0 0 

( \ - ^s) 


1 t 


(6.26) 


ji ~ 2 > i - 2, , (6.27 ) 

and are [M ] - orthonormal vectors. When a partial eigen 
solution is available, the derivatives of an eigen vector can 
be approximated as 





\ok 


(6.28) 


■where p < U. 


This approach has boon used by Kapoor to find the eigen 

vectors and dynamic dispLacomonts at modified designs . Recently, 

77 

Imani, Sandor and Kramor used this approach in thu deflection and 
stress analysis of high speed planar mechanisms V7ith elastic links. 


The quantities 


sS] 

ac. 

0 


and 


liM 

ax. 

0 


are obtained in the 


same way as explained in the pre-vious sections. In the present 
i7ork, using tho equation 6.28, the approximate eigen vectors are 
compatod at 1 and then tho corresponding values of are 

computed. The values thus computed are compared with those obtained 
by an exact oigen solution in Table 6.1 . Here again, it can be 
observed that tho two sets of values agree well -with ea,ch other. 



TAH>3 6 .1 COlIPAxilSIOIT OF RLSPOKSE QU/J^TITIES BY EXACT iJ® APEROXILIATE ISTHOD 
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FIG.6-1 COMPARISION OF FUNDAMENTAL NATURAL- 

FREQUENCY BY EXACT AND APPROXIMATE METHOD. 



CIIiPi'liiR 7 


ILLIIST-RATITE EXMiPLES 

TVie optimization protlein foirtnalated in the previoas clmpte 
ho-s been applied to the miniiaim weight design of two different types 
of machine tool structures, naLiely, the lathe bed and a horizontal 
knoe-type of milling machine structure . These machine tool 
structures were optimized in order to demonstra.te the feasibility 
and offectiveness of designing complex structures '.•■ith multiple 
behavior constraints. A computer program has been developed to 
o].vbimiz'e machine tool structures using finite element idealization. 
As stated earlier, the analysis routine incorporates the frame 
elements, triangle bending elements and triangular piano stress 
elements. These elements -are sufficiently general to idealize 
several types of machine tool structures. The static and dynaijio 
analysis programs as well as xhe optimization program are also 
written as general subroutines. Hence, the optimum design of any 
other type of machine tool structure can he accomplished by using 
the present comjutor program hy making little modifications. The 
computer piogram developed is described in Appendix B. 

The results of the present optimization study are presented 

in this chapter. 

7.1 LA'J™ BED BHSIGN 

Ex‘;.iriple 1 ; 

The first example considered is a Warren-type lathe bed 

all diamensions of this bed are taken 


shown in Fig. 7 *1 • over-i 



from the bed of lo,tb.e machine model lo. MB02, year , manufac- 
tured by Mysore Kirloskar Ltd., Harihar, India. This bed is 
optimized with a lower limit of 65O radians/ second on the funda- 
mental freq_uency and an upper limit of 0.00174 radians on the 
torsional deflection when a torque of 15,000 kg. cm. is applied at 
one v,.nd of the structure by keeping the other end fixed. This 
torque corresponds to a horse power of 4 iT ^ke turning is done 
at tho lowest spindle speed of 22 rpm. In this design, 4 design 
variables and 10 design constraints are considered. Tlie finite 
oloiaent modelling for this structure is shown in Fig. 7*2. The 
idealization of the structure consists of 26 node points, 46 
triangular plate bending elements and 52 frame elements. The 
number of degrees of freedom considered in the static and eigen- 
value analysis is 152. The band width of the stiffness matrix 
is 56. The guides on the main members arc taken to be rectangular 
in cross section(6.5 cm x 2.0 cm.). 

In the generation of the global stiffness and inass 
matrices, identically disposed finite elements of tho same size 
and shape are grouped together so as to reduce the computational 
effort involved. This resulted in the saving of computer time 
upto about P,Ofo in the generation of global stiffness and mass 
matrices. The design variables are 3 

1 

= thickness of the main members 
X2 = thickness of the lacing diagonals 
X = stiffener width on the main and lacing diagonals 
X = stiffener depth on the main and lacing diagonals. 
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The mterial of the atruoturo is taken as gre;r cast iron v/ith the 
following material properties s 

density (P ) =7.2 gms./cra^. 

Young's modulus (s) = 1.12 x 10^ Icg/cm^. 

Poisson's ratio (v) = O.25 

The geomotrical and hehavior constraints are as follows *. 
Geometrical constraints 



= 

1.0 

cm. - 

-^1 

§2 


V 

^“1 

' 2.0 

cm. 

% 

= 

1.0 

cm. - 

^4 


h - 

2.0 

cm. 



1.75 

cm. 

- X. 

% 

= 

3- 

■ 3.2 

cm. 

gy 

= 

1.0 

cm, • 

-^4 


= 

h- 

‘ 2.2 

cm m 


Behavior constraints ; 

g^ = d^ ^ 0.00174 radians 

®10 ~ radians/sccond - co^ 

where d^ is the average torsional deflection of one end of the 
hod in radians. The optimization remits are shown in Tahle 
The progress of the optimization path, showing the cumulative 
nuinhor of one -dimensional minimizations verses the weight of the 
structure, for this example is given in Pig. 7 *3* The least 
weight design has a weight of 47.93 hg. with a rcjducticn of 
45.5'/!> corapared to the starting design. 
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Among the side constraints, and arc at their 
lower hounds. Ilono of the behavior constraints are active in this 
exo,mple. The natural frequency gradually increased from 
71*1 radians/second at initial design to 755 radians/second at 
optimum design, whereas, <02 gradually decreased from 978 radians/ 
second at initial design to 949 radians/second at optimum design. 

The average computer tine taken for one one -dimensional minimization 
is 3«5 mihutes on IBi 570/l55 computer. The nutiibor of ono- 
diiaonsional minimizations is 10 and the computer time taken is 35 
minutes. 

Example 2 s 

In this example, the same lathe bed shown in Pig. 7»1 
is optimized by taking 5 design variables and 14 design constraints. 
The design variables X^ , X 2 , X^ and X^ arc sane as before. The 
vfidth of the lathe bed is included as the fifth design variable 
X^ in this example. A limitation is. also placed on the second 
natural frequency in this example. The initial design vector is 
chosen such that the behavior constraints arc nearly satisfied 
at X^ . This was done to see whether any weight could be reduced 
from this starting design vector. Tlie starting design corresponds 
to a weight of 87*49 ^^* geometrical and the behavior 

constraints in this example are as follows : 

Geometrical constraints j 
= 1 .0 cm. - X., 

= X.| - 2.0 cn. 
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~ 1.0 cn. — X2 

£4 == ^2 - 2.0 cm. 

§5 = 1 -75 cm. - 

Sg = “ 3*2 CEi. 

gy = 1.0 cn. - 

gg = X^ - 2 .2 cm. 

g^Q = 10.0 cm. - X^ 

g^^ = X^ - 18.0 cm. 

Behavior constraints ; 
g^2 ~ ~ radians 

g^j = 700 radians/s 000 nd - 

§14 = + 100 - Wg 

'Iho optimization results arc shown in 'lahle 7*2. The 
torsional deflection 'constraint is active in this example. The 
width of lathe hed (X^) at optimum design is 13*50 cn. and it 
corresponds to a lacing angle- of 55 •5*’ • ^he proposed least 
'vYoight design has a w'eight of 77*42 kg. with a reduction of 11. 45^^° 
coraparod to the starting design. The progress of the optimization 
path is shown in Fig. 7.4. The number of ons-dimonsional 
minimizations is 14 and the computer time taken is 55 minutes. 

7 .2 HORIZOIJTAL KHEB-TYFE MILLING MiCHENE STEBCTUSIL DESIGN 
Example 3 (a.) : 

The horizontal knee-type milling machine structure, 
shovm in Fig. 7.5, is considered for optimization in this example. 
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The finite element modelling for this stractnre is shown in Pig,7»^» 
In this example, 6 design variables and 16 design constraints are 
considered. The idealization of the structiir<. consists of 50 node 
points, 50 triangular plate bending elements ^,..d 18 frame elements. 
The number of elastic degrees of freedom* is 5- a-nd the band, width 
of the stiffness matrix is 42. In this problem, about 20fo of the 
computer time was saved in generating the global stiffness and 
mass matrices by grouping finite elements of the saiae size and 
shape having the same transformation matrix from local to global 
coordinate system. 


The arbor diameter is taken as 4 arbor support 

is idealized as two frame elements of rc-ctangular cross-section 
(10 cm. X 7 cm.). The thickness of the dove-tail on front face 
of the column is taken as J.Scm.and it is added to the thickness 
of tho column (X^) on the front face. The thickness of the dove- 
tail on tho bottom side of the ovor-airm is taken as 2,5 cm. This 
thickness is added to the thickness of the over-ara (X^) on the 
bottom side. 


The design variables are s 






column depth at the bottom 
thickness of the over-arm 
width of tho machine 
thickness of the column and table 
column depth at the top 

square cross-sectional diamonsions of the ribs on 
the over-arm and its joint wiuh the column. 
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The geometrical and the behavior constraints aro as follo¥;s ; 
Geometrical constraints' ; 



= 

- 54*0 cm. 

S2 


35*0 cm. - X^ 

§3 

= 

0.8 cm. - X^ 



X2 - 4-.0 cm. 



24.0 cm. - 

^6 


X^ - 35*0 cmv 


= 

0,8 cm. - X^ 

% 

= 

X^ - 4,0 cm. 

gg 


^ 5 -"l 

®10 

= 

30.0 cm. - Xj. 

5 

®11 


-^6 

^12 

= 

Xg - 6.0 cm. 

Behavior constraints s 


= 

d - 0.009 cm. 

G 

®14 

= 

IGmIjjI - 0.000003 cm./kg. 

gl 5 


850.0 radians/ second - 

®16 


+ 100 - “2 


whore^ 



laaxinain cutter centre deflection, and 




minimm negative inphaso cross receptanco of 
cutter centre relative to table. 
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The material of the columj table 5 over-^arn and arbor 
support is taken as grey cast iron with the samo material properties 
given in the case of ^.’^i^rren typo lathe hod. The material for arbor 
is assumed as \/r ought steel with the following material properties ; 

donsity (p) = 0.78 gas./cmb 

g ? 

Youiig's modulus (E) = 2.1 x 10 kg./cm . 

Poisson's ratio (v) = 0.3 

The optimization results are tatulated in Table 7*5* The 

starting design correspemds to a weight of 74O kg. The proposed 

optimum design corresponds to a weight of 3OO.78 kg, with a 

reduction of 45 •5/'° weight compared to the initial design. The 

active behavior constraints are static deflection (d = 0.006932 cm,) 

c 

and the difference of and ” ^2 ” radians/second). 

The constraints on and approached criticality as the 

Optimization process progressed. Among the geometrical constraints, 
Xo = -987 cm. and X = O.905 cm. arc near to thoix respective lower 
bounds. The average computer time required for one one-dimensional 
viinimization is 8 minutes on lai 370/l55 system. The progress of 
the optimization path as a plot of 'f and 'p' functions verses 
the cumulative number of ono-dimensional minimizations is shown 
in Fig. 7 '7 • 

Example ^{'b) t 

In order to see whether the optimum design obtained in 
example 5(a) corresponds to a local minimn or the absolute minimuia 
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in the design space, the satie example has been renin with a 
different starting design vector. The plot of 'f and d- 
fanction as the optimization progressed is shown in Fig. 7*8. 

It can ho observed that the plot is similar to the one shown 
in Fig. 7*7* The optimization results for the oxaiiplo are shown 
in Table 7*4* The optimum, design variables in the two cases 
a-gree well with each other except for small differences that might 
have occurred due to some round of errors a.nd numerical 
instabilities in the optimization process. Although, merely on 
the basis of two trial starting designs, it is hard tc say that 
the minimum obtained is the absolute minimum over the design space | 
finding the similar least weight design by starting from two 
different initial designs is atleast a pointer in that direction. 

7 .3 SENSITIVITY ANALYSIS 

In practice , a designer \/culd be interested in knoiving how 
the resp'mse quantities vary with a change in the design parameters. 
This type of sensitivity .analysis will help the designer in ir^nipu- 
lating the design variables to suit sene specific requirements- 
Jhrthor, in some Cfases, the results obtained from the optimization 
procedure may have to be rounded off tc .the nearest practical values 
of the design va.riablos. Hence a sensitivity analysis -'f the 
response qua,ntities d^ , a.nd with respect to the various 

design variables is corductod. In this analysis, the reference 
design is taken same as the starting design point of example 5* 
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The results cf the sensitivity analysis arc shewn in Pigs. 7 *9 
to 7«11 • Pron these figures, it can bo seen that all the 
response quantities are most sensitive with respect to the 
vfidth of the machine. The responses d^, and can bo 

observed to Ki least affected by a change in the oross-section 
of the ri.bs on the ovcr-aimi , The natural frequency ha.s 

been found to be least affected by a change in the depth of 
the column at the ba-se. It cp.n also be seen that the fundaraental 
frequency varies nearly linearly v;ith respect to each of the 
desi^^pn variables. 
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•AGTrVE CONSTEAIHTS 





TABIB 7 ^3 OPTIMIZATION RESULTS OP HORIZONT/iL MILLING MACHBiB (EmiPLl 3 a) 
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CHAirSE 8 


EJaiABILITY ISiAlYSIS OP MACHIHE TOOI SIHJCTO'EES 

mie conventional approach, in design practice, which considers 
loads, resistances and design variables as deteiministic quantities,, 
may be compared to a kind of worst case analysis. In this ^proach, 
as an exaiiiple, when designing a structure for strength, the maximum 
of loading and the minimum of strength are treated not only as 
representative of design situations, but also of simultaneous 
occurrence. Ihe design of a structure for an extremely large load 
(of rare occurrence) and for an extremely low value of strength 
(of rare incidence) will be very uneconomical. In such situations, 
the loads and strengths have to be treated as random variables and 
the probability of simultaneous occurrence of extreme values of 
load and strength becomes very important. 

In probability (reliability) based design, usually a 

design is considered to be safe and adequate if the probability 

of failure of the structure is less than a specified small quantity, 
“6 

say, 10 . Ihus the basic step involved in reliability baaed 

design is to analyse the structure for reliability. In this chapter, 
a method of analysing the reliability of machine tool structures 
is developed. More specifically, the reliability analysis of 
horizontal milling machines in various failure modes is considered 
for illustrating the method. 
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8.1 cnoics 0^ EAlfflOM TAHIiBLES 

The operating conditions, material properties and modal 
daraping factors vary substantially in machine tool structures- 
In the present work, the r.;liability of the milling machine 
structure considered in example 3 of chapter 7 is evaluated at 
the initial and optimum designs. !Hie reliability of the structure 
against the various response quantities of failure modes is 
found by taking the table height (h), the distance of tne cutter 
centre from arbor support (c), the damping factor (?), the Young’s 
modulus of the material of the structure (e), and the load on 
the cutter and the table (p) as random variables. The response 
quantities considered are ; ( 1 ) the maximum static deflection of 
tne cutter centre in any direction, (2) the first natural frequency 
of vibration, and (3) the minimum negative inphase cross- 
receptance of the cutter centre relative to the table. 

The Location of the table and the cutter centre depend 
on the size of the work piece. Since the size of the work piece 
will vary for different jobs, h and c are taken as random variables. 
Erom the present knowledge, the modal damping factors of structures, 
particuleirly in machine tool structures where joints are involved, 
cannot be estimated precisely. Therefore modal damping factors 

are taken as random variaoles. In actual practice the material 
properties vary and hence the Young's modulus (e) is considered 
as a random variable in this work.^ linally the magnitude of the 
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Ilhe probability of failure is given by 


= P(e < i) = /“ . f(jl) . d£ 

—00 


1 - / 


( 8 . 1 ) 


v/here f„(x) and represent liie probability density and 

jL jl 

distribution functions respectively. 


If several loads act simultaneously on the structural 
system as shown in fig. 8.2(a), the failure probability is given 

fey 

OD P 

= 1 - / [ n (r)] fj^(r) . dr (8^2) 

d=1 3 

If a single member is subjected to several load conditions as 
snovfti in Fig. 8.2(b), the probability of failure can be determined 
from the relation 

00 1 

if = 1 - / { n [ 1 - F (\) 3 } f^(£) . d£ (8.3) 

_oo i=:1 

where £ . is the force induced in i'*'^ member due to the load Ji- . 

1 

Since this is a vreakest-link chain, the term in braces in 
Equation f'8.3 represents the probability of survival of the chain 
and is based on all links surviving under the load I = £ . Hiis 
term is evaluated from the products of probabilities of individual 
links surviving under the load L = £ . Equation 8.3 is often 
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approximated as 


/5. 1 - n (l - P. ) 

i=1 i 


(8.4) 


.th 


where denotes the probahility of failure of i ' link, iinally 
^i 

the multicomponent, multiload system shown in Pig. 8.2(c) is given 
by 


If = 1.- / / ••• / 


h ^2 ^ ^=1 


^^2^ *•* h . d£2 ••• d£ (8.5) 


th 

where £. is the total force induced in i member due to the 
am 

lo.ods jl . Since the computation of the exact probability 

12 p 

of failure is a complex prob>.iDilistic problem, Equation 8,4 is often 
approximated as 


q P 

P 1 - E n (1 - 1 ) 

i=l D=1 


( 8 . 6 ) 


th 

where P^ denotes the probability of failure of i member ynder 


.th 
j load 




8.3 CUStE'UTATIOS OP RELlABIlITr IK A PARIIOJLAE PAlIUES MODE 

The reliability (R^) of a system is taken as one minus the 
probability of failure (P^). The method of computing the reliability 
of a system in any particular failure mode is considered in this section. 
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If R is the resistance and I is the load acting in the specified 
failure mode the reliability of the system can be analysed as a 
single member-single load problem. For simplicity the resistaiice 
and the load are assumed to be normally distributed so that 




and 




1 . r 1 /r - En 2 i 

— exp [- - {—^) i 

R 


(8-7) 


( 8 . 8 ) 


where L and E represent the mean values and cr arid cr th$ 
standard diviations of I and E respectively. althou^ Equation 8 . 1 
is applicable, the following simpler procedure is used to find the 
reliability of the system, in this case- S'/ defining a new random 
variable C as 


5 = R - i , 


(8.9) 


the reliability of the system can be expressed as 


E = P(5 1 O) = /“ f^ (n) dT) 


( 8 . 10 ) 


where f^(h) is the density function of E given by 


1 


f ('I) = 

^ /2'it 




— exp L p 

a ^ f. 

o_ s 


( 8 . 11 ) 


where 5 is tne mean value and is the standard deviation of ^ 
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If H and L are independent, 
given by 


the expressions for 


g and 0 ^ 


are 


5 = E - L 


( 8 - 12 ) 


and 


2 2 
a_ =/0_ + a_ 


(8.13) 


Equation 8.10 and 8.11 give 


1 [ r 1 /h - 

E = — j exp [- - (—5—) ] dn 

° /2ir 0 0 ^ ? 

5 


(8.14) 


By defining a standard normal variate Z as 


Z - n - ? 

'■“s’ 


Equation 8.14 can be rewritten as 


(8.15) 


Rq = "i: r i 2^) <12 

Z=z 


( 8 . 16 ) 


where the lower limit of integration (z) is given oy 


z 



E - I 


2 

I 


] 


(8.17) 


Once the value of z is calculated, the corresponding reliability 

can be determined from equation 8.16. This value can be obtained 

78 

more readily from the standard normal tables . 
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In order to find the reliability of a structural system. 

using Equations 8.17 and 8. 16> the mean values and standard 

deviations of the generalized load and the generalized 

resistance (h) must be known. Since E and E generally depend 

on several other randan design parameters, one has to determine 

E, E, 0^ and in teims of the means and standard deviations 
Xj r 

of the random design parameters. In general, if Y is a nonlinear 
function of several random variables x. , x , ...,x , the approximate 

12 . S 

values of the mean and the variance of Y can be found by linearizing 
Y about the mean values of x^, using Taylor’s’ series 

expansion. The expressions of Y and are given by 


and 


^ *2 V 


S 

Oy % ^ I 
1=1 


[il-l - _ 

3x^1(x^,X2, 



(8.16) 

o^ (8.19) 

X. 

i 


v/here the random variables x. jX , •..,x are assumed to have zero 

1 z s 

correlation. 

Since the reliabilities against the response quantities 

d , 0 )^ and are to be found, the partial derivatives of d , 

0 1 Mil" 0 

( 0 .J and with respect to the random design parameters h,c,C»E 

and P (evaluated at the mean values of the design variables) are 
required. These partial derivatives are found by using a finite 
difference scheme in this work. 
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8.4 ItJi£EIRICiiL EESDIIS 

21he reliability of the horizontal milliiig machine structure 
considered in example 3 of the previous ch£^)ter is analyzed in 
each of the failure modes. !I3ae partial durivatives of the 
cutter centre deflection, fundamental natural frequency and the 
minimum negative cross-recej.tance of cutter centre relative to 
the table have been computed by finite difference scheme by 
talcing a 5^ variation in the design parameters. 

!Ihe reliabilities of the machine tool structure are 

found at the initial and optimum designs of exnriple 3 of chapter 7. 

The values of the partial derivatives of the response quantities 

with respect to the design parameters are given in Table 8.1, 

The reliability results obtained by using a coefficient of 

variation (v ) of ©.05 for all the random design parameters 
i 

are shown in Table 8.2. Iroa these results, it can be seen 
tnat the reliabilities are lower at the optimum design compared 
to tnose at the initial design. This is to be expected in this 
case since the optimum design has been obtained according to 
deterministic formulation and most of the behaviour constraints 
(involving the response quantities) terri to reach criticality at 
the optimum point. If the optimum design is obtained according 
to probabilistic formlation (using Equations 2.2 of chapter 2), 


the reliabilities corresponding to the optimuii design are expected 
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to be higher than those corresponding to the initial design. 

In many practical problems, one would be interested not 

only in the reliability of a structure for a particul.or set of 

data, but also in knowing how the reliability changes with a 

variation in the standard deviations of the various design 

parameters, ' Hence the variation of reliabilities of the milling 

machine structure against d^, and various coefficients 

of variation of the design parameters (v ) are shown in iigs.8,3, 

^i 

8.4- and 8.5 respectively. H'ly, these figures, it can be seen 
that the reliability in the deflection failure mode is most 
sensitive with respect to and least with respect to at 
both initial and optimum designs. Ihe reliability against the 
fundamental frequency is most sensitive to a change in Vj, and 
least to a change in and Similarly, the reliability 

against the response quantity can be observed to be most 

sensitive to and least to Vp. 

If the machine tool structure is considered as a weakest- 
lihk chain, each link representing one particular failure mode, 
the overall reliability of the system can be calculated from 
Equation 8,4. This gives the overall system reliability as 
(r„) = 0.99895646 88 at the initial design and 

(^o)overall ~ 0.4816022724 at Ihe optimum design. 
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8.2 EESULTS OB' EELI.JBILITy .JJ.ILYSIS mTE RESPECT TO THE VARIOUS RESPONSE QUANTITIE 
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Note : Coefficients of v..riatioa of all the design variables = 0.05. 
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CEAPTER 9 


Discussion AID GOICIUSIOIS 

The results of ite example problems presented in chapter 7 
demonstrate the feasibility of automated optinaim design of machine 
tool structures with static, dynamic and regenerative chatter 
stability constraints.. The numerical results of chapter 8 indicate 
the feasibility of reliability analysis of complex machine tool 
structures. In this chapter, the results given in the previous 
two chapters are analysed and the conclusions drawn therefrom are 
enumerated. In designing these complex machine tool structures 
with multiple behavior constraints sevsrsl analyfeisi programs have, 
been developed and incorporated into the optimization routine and 
therefore it is difficult to summarize all the findings. lEie 
conclusions drawn from the present study a 3 ?e grouped in section 9.1. 
Rinally suggestions to extend the scope of the design problems are 
given in section 9.2. 

9.1 SOMMiiRY .tiED COICLUSIOIS 

(l) In the finite element idealisation of ¥arren type latho 
bed and horizontal milling machine structures, triangular 
plate bending elements and frame elements have been fomd 
to be simple and efficient. A moderate number of finite 
elements have yielded reasonably accurate results i She grouping 



145 


of finite elements of the same sise, haying the same 
transformation matrix from elemental to global coordir.ate 
system, has resulted in a saying of about 20 to 8C^ of the 
computer time in- generating the global stiffness and ma s 
matrices. In the reanalysis cycle of optimisation procedure, 
the generation of stiffness and mass matrices of the 
structure is required a number of times and hence the grouping 
of such elements has resulted in saying of considerable computer 
time. 

(2) Ihe cholesky decomposition of symmetric banded matrices has 
efficiently been used in solving equilibrium equations and 
in obtaining a partial solution of eigenyalue problem by 
Ralei^-Ritz sub-space iteration algorithm. This algorithm 
solves the eigenyalue problem directly without a transformation 
to the standard form. This method is the most efficient one 

for solving a partial eigen solution of large structural systems. 

(3) Ihe receptances of the cutter centre relative to the table 
of horizontal milling machine have been obtained by usi 'S 
modal coordinates assuming the damping matrix as a linear 
combination of stiffness and mass matrices. 

(4) In this work, the geometrical constraints are linear, vdiereas 
the behavior constraints and the objective function are non- 
linear. Hence the constrained optimum design problem is 
cast as a nonlinear mathematical programming problem. The 
computational experience shows that the approximate methods used 
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in the present work for evaluating the gradient and t;- 3 
slope of 4 -function have been quite efficient and reliable 
without involving any significant errors. 33ie progress of 
optimization has been quite smooth without any undue i xaber 
of optimization steps. 

(5) In computing the values of for regenerative chat-&er 

stability in designing milling machine structures, all the 
possible common cutting conditions have been considered. Ihr 
an unit amplitude of harmonic forcing, the negative inphase 
cross receptances have been computed at an interval of 

10 radians/second beginning from the first natural frequency 
to fourth natural frequency of vibration and from these, the 

mi n.i.:',iuti value (&. has been found out. 

^ MIIT 

(6) Ihe coupling constants in plain milling operation are 
required in the design of horizaatal milling machine structures 
for regenerative chatter stability. Ihe coupling constants 
are obtained, for most common range of cutting conditions, 
both by the application of basic theory of metal cutting and 
by using a procedure based on an empirical formula suggested 
by Vulf. l!he values obtained from both these methods agree 
well with each other. The limiting value of the coupling 
constant (r^^,,) has been selected from this analysis. 

(7) Prom the optimum results of lathe bed (example 1), it has 
been found that the thickness of both the main members and 
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lacing diagonals decrease as the optimiaation progressed. 

At the optinum, the thicknesses are found to he the same and 
are at their lower botuads. iHiis indicates that the thin 
structures are preferable for this type of lathe beds. 2iie 
widths and thicknesses of flanges (stiffners) on main maabers 
and lacing diagonals also decreased as the optimization 
progressed. In example 2, the width of the lathe bed is also 
included as a design variable acd the behavior constraints 
are taken as nearly active at the initial design. !Ilhe results 
show that the width of the lathe bed increased with a corres- 
ponding value of the lacing angle as 53.5°., Ihis value of the 
lacing angle agrees very well with the results obtained by 
Badauri, Moshin and Ihornley who analysed and tested various 
■Warren beams under static torsional loads. It may also be 
mentioned that the thicknesses of both main members and 
lacing diagonals decreased approximately at the same rate 
and at the optimum, their values were same. Iherefore, it 
c.-jn be concluded that it is preferable to have equal values 
of thickness for the rnaiij. matbers and lacing diagonals. In 
designing this type of structures, the machine tool structural 
designer can adopt this as a general guide line. 

(s) Ihe optimization results of the milling machine structure 

(exaaples 3o, and 3b) show that the thicknesses of the over-arm, 
the coluiiin and the table decreased as the optimization progressed. 
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Eyen, in horizontal milling machin. structures, it can be 
concluded liiat thin structures are preferable. !I5ie cross 
sectional areas of the square ribs on the oTer-aim and its 
joint with the colujan also decreased as the optimization 
progressed. For Ihe same width a.id depth of cut in plain 
milling, the static forces on the machine can be reduced by 
decreasing the feed. Hence the design becoaes safer from 
the static rigidity point of view. However, the sane can’t 
be said about the chatter stability of the machine. Before 
the rated maximum horse power capacity of the machine is 
reached, the regenerative chatter may develop depend ong upon 
the orientation of the cutting forces and the value of the 
corresponding coupling constant. Since chatter stability 
detenJines the metal removal rate, it becanes essential that 
the value of the coupling constant be taken for the worst 
case in designing milling ’machine structures for chatter 
stability. 

( 9 ) At initial and optima i designs, the reliability analysis of 
milling machine structures in various failure modes by taking 
the average height of the table (h), position of the cutter 
centre (c), modal duiaping factors (?^), Young's modulus (e) 
and the loads (b) on the machine as probabilistic quantities ^ 
Ihe results shovf that at optimum design, the reliabilities 
are veiy much lower compared to those at the initial design.. 
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Ihis is expected,, since the behavior quantities are far 
fron the allov'able lioits at the initial design, iccoxding 
to the results obtained, at initial and optiaua designs, 
the reliabilitj against deflection is more sensitive vvith 
respect to a and a_,. Ihe reliability against deflection 
wiih respect to and a is not as sensitive as with 

xl C 

respect to and c^,. Eie reliabilities of the structure, 
at initial and optinum designs, against a),j are ncre sensitive 
with respect to Cj, and With respect to the variation 

of other probabilistic quantities (a , a_ and 0 ), the 

c ^i ^ 

reliabilities do not change. Ihe reliabilities of the 
structure at optiaun design against that the 

sensitivities vd.ll increase with respect to the coefficients 
of variation of the probabilistic quantities 5, c, 3 
and h in increasing order. At the initial design, the 
reliability against reaained almost constant at 1. 

9.2 EECOmiMDATIOIS BOB BUREHER TOHK 

( 1 ) Im obvious extension of the present work is to apply the 
procedure of amlysis and autotiated opti-iun desiga to 
the complete lathe machine structure including the gear 
box casing, the spindle and its supports. The location 
and the number of spindle supports and their flexibilities 


can be included as design variables. 
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(2) In a adlling machine stiucture, the overall diaensious liice 
the depth of the column at the base, the depth of the column 
at the top and the width of the machine are very important. 

It IS difficult to summarize the general relationship between 
these major dimensions from the design example considered in 
the present woik. !Eherefore, to serve as a general guide for 
a machine tool stiuctural designer, it is worthwhile to optimize 
a number of horizontal milling machine structures for various 
borse power capacities, sizes of woik piece and limits on 
behavior constraints and establish a general relationship 
between these major dimensions. 

( 3 ) In horizontal milling machines the over-arm forms an important 
part and, therefore, for a given colunn of the milling machine 
stricture, the over-ar... can be further optimized by taking 
the following qu;antities as design variables : 

i) number of ribs on the stiuctuie 

ii) cross-sectional shape on the ribs 

ill) taper of the over-arm on the sides 

iv) taper of the over-arm on the top 

v) depth of the over— arm 

Vi) widtn of the over-.urm. 

rne analysis c.'-'n be made either by taking the joint of the 
ovo3>-^r-.i with the column as a flexible support or by idealising 
the column as a whole by frame elements. 
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(4) Hie design of -machine tool structures for other criteria 
like maximization of static rigidity, fundamental natural 
frequency and regenerative chatter stability can also be 
studied, 

( 5 ) Ihe reliability analysis procedure presented in this work 
can be extended to cases where the design parameters follow 
distributions other than normal distribution. 

(6) The operating conditions, material and structural properties 
(damping factors are influenced by the joints of the structure) 
vary substantially in machine tool structures. Instead of 
optimizing machine tool structures by taking the most 
unfavourable cutting conditions and sizes of work piece, the 
reliability based optimum design of machine tool structures 

can be attempted to get a more realistic design by considering 
the operating conditions, material and structural properties 
as probabilistic quantities. 
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APEETOIX A 


ELHmTT STIPHffilSS AUD MASS 1ATEIC13 


A.1 STimiESS MD IIASS MTRICES of a TPJAtTGULAR PLIi.TE 

EBKDHTG ELEWT 

In deriving the elemental stiffaess matrix [K] for a 

triangular plate "bending element (shown in Fig. A. 1), when "both 

in nlane and "bending stif fn'’-'?8es are considered, equations 5 ‘If 

5.2 and 3*3 of chapter 5 used. From equations 3*'I ond 3.2, 

the in-plane stiffness matrix ( lO of the element will result in 

P 

a 6 X 6 matrix in loccA coordinate syston. The me.trix . 

^ Jr 

can "be written as follows; 


!>^lp 

(6x£) 


(2x2) 

(2x2) 

(2x2) 


^12^ p 
(2x2) 

^ ^22^ p 

(2x2) 

f'32^P 

(2x2) 


t^13ip 

(2x2) 

^ ^23^ p 
(2x2) 

^^35’p 

(2x2) 


(A.1) 


J 


i’j;}_..pl 3 ,ne analysis, the relationship "between the node.! dis_ 
placements and forces of an element can "be written as ■ 
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whore u.^ and indicate the oomponents of displa-cement of node 

i(i= 1 , 2 , 5 ) parallel to the local x-and y-*axes respectively. 

Similarly P and P denote the^omponents of foi»e at node 
i ^i 

i(i= 1 , 2 , 5 ) parallel to the x-and y-axes respectively. Hotice 
that the suffix ' p’ has been used to indieate the in- plane 
stiffness matrix in equation 

Similarly a relation "between the forces and displacomaits 
corresponding to the "bending of the plate (obtained from equation 
■5.5 of chapter 3) can be TOitten as ; 








(A.3) 
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■vrhere [K]^ is the bending stiffiiess matrix in the local 

coordinate system. Here and Indicate the components 

of displacement and force parallel to z axis at node i, and 

M and M represent the generalized forces corresponding to 
i ^i 

&•_ and0„ at node i respectively (i s= 1,2,5) the suffix 'h' 
^i 

has been used to denote the bending stiffness matrix. 

The order of the bending stiffness matrix, , is 

9x9 a^d it <:an be written in the following form : 



In the analysis of three dimensional structures, both in-plane 
and bending stiffnesses have to be combined. In combining these 
stiffnesses the following points are talcen into account; 

i) For small disple-coments, the in-plane and bending 
stiffnesses are uncoupled, 

and ii) the in-plane rotation 6^ (rotation about the local 
z- axis) is not necessary for a single element. 



However, 6^ and its conjugate force have to he 
considered in the analysis hy including appropriate 
numher of zeroes to obtain the elemental stiffhess 
matrix (for the purpose of assembling several elements) 
Therefore to obtain the elemental stiffness matrix 
the in-plane and the bending stiffhesses are combined 
as shown below ; 


_ . - _ , I II _ . 4 

0 ’ 0 

I 

0 0 




0 


0 0<0 


[SoJ V ' Oi 0 oi 0 [s^,] 


25 P } 

lO 0 

0 o'! 


0 ot 


0 0 0 0 0 ’ 0 : 0 0 iO 0' 

I ? 


0 * 0 
t 

0 ^ 0 
I 0 

* 0 
i 

I 0 


t 

0 I 0 


52 p 

















165 


A. 2 nMssmmmioi to global coohduiite systei! 


The stiffness matrix derived above (equation A. 5 ) 
is ¥vith reference to the local coordinate system shown in 
Fig. A.1 . In ■file analysis of complex structures having 
finite elements with different orientations, it becomes 
necessary to transform the local stiffness matrix to a common 
set of coordinates Imown as the global coordinate system* 

This can be achieved with the help of a transformation 
matrix, [A] . The global stifftiess matrix, , of an 

element can be written as s 

T 

Ek]q = [A] [K] [A] (A. 6) 


1 

where [A] represents the transpose of [A] 
The transformation matrix is given by 


fr 


1 


[A] 

(in X 1R) 


'"I ^ , 

[!■] 

Cl 


whore [I] is given by 


m 

(6x6) 




n. 


^3 

0 0 0 


0 

0 


0 

0 


0 

0 


J 


m^ n^ 0 0 0 


0 0 0 I 


1 1 1 


”2 ^ 


C ”3 




(A. 7) 


(A.8) 



wh6re , HLj and n^ represent the direction cosines of the 
local x-direction 'with respect to the X, Y stnd Z direction 
of the global Coordinate system respectively. Similarly the 
suffixes 't* and '5» represent the direction cosines of the local 
y and z directions. 

The elemental (consistent) mass matrix can he derived 
in a similar manner. However, in the present work, the lumped 
mass matrix, hy neglecting the rotary inertia properties, is used 
for obtaining the eigen solution. This simplifying approximation 
effected considerable savings in compiter storage ani time at the 
exT>ense of insignificant loss of accuracy in the eigen solution: 

A 5 HTiEi'IENTAL STIPFIIESS AHD J/h\3S MilTEICES OP A PElftME ELTMMT 

The elemental stiffness and mass matrices of a frame 
element given in reference 40 have been used in the present work. 
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A.1 IN-PLANE AND BENDING DISPLACEMENT IN 
LOCAL COORDINATE SYSTEM 


APPEKDIX B 


DESCRIPTION OF THE COMHJTER mOGRAli 

The main computer program for the optimization of 
machine tool structures is written in PORTRMT lY language 
it consists of 13 subroutines. Initially the computer program 
is perfected on IHI 7044 (s-t Indian Institute of Technology, 
Kanpur) and the final optimization results are obtained on 
IB5J 570/155 computer at Indian Institute 0 f Technology, Madras. 

B.1 PURPOSE OP 'THE SUBROUTINES 

(1) GSH ; To generate the stiffness and mass matrices 

for a triangular plate bending element and 
frame element in a banded form. 

(2) CBATU)? To obtain the cholesky decomposition of 

symmetric banded matrices. This routine 
stores only the upper triangular hand of the 
decomposed matrix and the diagonal of the 
matrix is stored in the first column, 
f 5) ZOLP ", To solve the equilibrium equations from 

the upper tria.ngalar band of the decomposed 


stiffness matrix. 
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(4) MSIC ; 


( 5) EIGEN ? 

( 6 ) 

(7) liilTIULs 
(=») BiOBULi 

(9) RESP 

( 1 0) ETN 2 

( 1 1 ) PTNZ 2 2 


(12) OBIiEVA; 

( 1 3 ) GTIMN f, 


(14) SGr..Q,UAs 


To obtain the partial eigen solntion of a 
structure by Rayleigh-Eit 2 s sub-space 
iteration technique. It calls routine 
EIGEN for solving the generalized Ritz 
problem , 

To compute the eigen values and eigen vectors 
of the generalized Ritz problem by power 
method. 

To find the inverse of a real matrix. 

To multiply two real matrices. 

To multiply two real banded matrices. 

To find the frequency response by the modal 
superposition technique. 

To compute the ^ - function by exact method . 

To compute the 0 -function by approximate 
method (making use of the expressions for the 
rates of changes of eigen values and eigen 
vectors) . 

To evaluate the objective function. 

To calculate the gradient of 0-function by 
evaluating the response quantities by 
approximate method . 

To evaluate scalar quantities for obtaining 
response quantities by rates of changes of 
eigen values and eigen vectors. 
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(15) STEESS s To oompate-ithe stress in different finite 
elements. 

(*^6) MTN 'TM4 s To implement the Davidon-Pletcher-Pos-ell 
varis-hle metric method of -unconstrained 
minimization. It calls the one dimensional 
minimization routine STEPX2 

(•tj) STEIX2 t It implements one -dimensional search method 
hy cuhic interpolation. 

T 

S'^jOPEZ j To evaluate the slops of function (S 



